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Abstract 



■^j- \ We introduce a new representation of generalized parton distributions and generalized distri- 
bution amplitudes that is based on the partial wave decomposition with respect to the complex 
^ . collinear conformal spin. This decomposition leads us to a versatile parameterization of these non- 
■ perturbative functions in terms of conformal moments, which are measurable for integer value on 
the lattice. This new representation has several advantages: basic properties and crossing rela- 
p tions are automatically implemented, a rather flexible parameterization is possible, the numerical 
Q_i! treatment of evolution is simple and analytic approximation of scattering amplitudes can be given. 

We demonstrate this for simple examples. In particular, phenomenological considerations indicate 
>► . that the t-dependence of Mellin moments is governed by Regge trajectories. The new representa- 

•rH , 

^ ■ tion is vital to push the analysis of deeply virtual Compton scattering to next-to-next-to-leading 
c3 ■ order. 
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1 Introduction 



Generalized parton distributions (GPDs) [U 121 El IH E] and their analog, obtained by crossing, 
generalized distribution amplitudes (GDAs) j2JIHl are non-perturbative functions that are accessi- 
ble in certain hard exclusive processes such as the hard electroproduction of photons and mesons 
off a nucleon or nucleus or hadron pair production by two photon fusion in e + e~ colliders. These 
functions are related to parton densities, form factors, and distribution amplitudes but contain 
additional non-perturbative information about the internal structure of hadrons and nuclei. Some 
of this information can even not be obtained in any other way than through a global fit of GPDs. 
This has been widely realized for the first time in connection with the proton spin puzzle. Here 
the second moment of a certain combination of GPDs provides the orbital angular momentum 
fraction carried by quarks of a given flavour |3] . Moreover, GPDs simultaneously possess a longi- 
tudinal and transverse momentum dependence and so it has been pointed out that they encode a 
three dimensional femto-holographic picture of the probed hadron or nucleus [Zj . Indeed, it could 
be shown that a partonic density interpretation holds in the infinite momentum frame as long as 
the longitudinal momentum fraction in the t-channel is vanishing (HI E] , see also Refs. JOl \H\ ■ 
More precisely, in the impact parameter space GPDs are interpreted as the probability to find a 
parton species % with momentum fraction x at a relative distance b± from the proton center. Even 
an interpretation of the three dimensional Fourier transform of GPDs in the rest frame has been 
suggested within the concept of phase space (Wigner) distributions jTSHTBj- For further details 
we refer to the comprehensive reviews in Ref . [Til ITS] . 

At present generalized parton distributions are one of the main topics of collider and fixed target 
experiments at DESY and JLAB. Further experiments are planed or proposed for COMPASS and 
ERIC. Also it should be mentioned that information on generalized distribution amplitudes has, 
e.g., been extracted from LEP data. Unfortunately, the wealth on information encoded in GPDs 
and GDAs goes along with their functional complexity. For instance, GPDs depend on both 
the momentum fractions in the s— and t— channel, x and 77, the momentum transfer squared t, 
the resolution scale Q, and the quantum numbers of the target and the probed parton. This 
multitude of functional dependencies is, however, very strongly constrained by their relation to 
parton densities, form factors, and distribution amplitudes and even more so by crossing relations, 
positivity bounds and Lorentz invariance in general. The latter implies in particular that the Mellin 
x-moments of GPDs must be polynomials of given order in the skewness parameter 77. 

Unfortunately, there is still another complication. Typically experimental observables allow 
only to determine convolutions including GPDs or GDAs and a formal deconvolution can prac- 
tically not be done for most of the processes 1 . To determine GPDs or GDAs from experimental 
data one therefore needs ansatze for them which involve only a minimal number of parameters. 

1 This is actually only possible if the hadron is probed with two virtual photons and the virtuality of both 
photons can be independently varied, which is experimentally an extremely challenging task [TBI fTTl ITS] . 
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Although, on the theoretical side GPDs have been intensively studied in the last few years, only a 
few of such parameterizations have been proposed and used in phenomenology. Perhaps the most 
popular parameterization is based on the ansatz suggested by Radyushkin, in which by construc- 
tion the relation to parton densities, form factors, and polynomiality is assured. However, this 
advantages arise from the simplicity of this ansatz which also might implement a certain rigidity. 
Especially, it is widely used in combination with a factorized t-dependence although the latter is 
known to be wrong. This factorization does not respect the disappearance of the t-dependence for 
x — > 1 120] an d is also basically ruled out by lattice results [2U 122 122] • How far the employed 
versions of this ansatz are suited for the kinematics accessible in present experiments remains an 
open question. Obviously, it is highly desired to have a versatile parameterization of GPDs and 
GDAs, which respects all of their formal properties. 

The main idea guiding the search for a more appropriate parameterization of GPDs and GDAs 
is that the relevant kinematic variables should be separated in this new representation. Let us 
remind how helpful the representation of parton densities by Mellin moments has proven to be for 
the analysis of hard inclusive processes. Mellin moments are given by the analytic continuation 
of the forward expectation values of leading twist-two operators with given spin J. The main 
advantage of the Mellin space is that operators with different spin J do not mix under evolution 
and so the solution of the evolution equations is trivial. In the case of GPDs and GDAs lead- 
ing twist-two operators can contain total derivatives and so the operator basis has to be chosen 
differently, in such a way that the operators again do not mix under evolution. The appropriate 
operator basis is given in terms of collinear conformal operators that are labelled by the (collinear) 
conformal spin and the normal spin of the operator. The former quantum number characterize 
the irreducible multiplets or conformal towers of the collinear conformal algebra while the latter 
denotes the members of a given multiplet. Group theoretically we are dealing with the representa- 
tion of the so-called collinear conformal algebra so(2, 1) which is a subalgebra of the full conformal 
algebra so(4, 2). Let us remark that except for the trace anomaly, which is proportional to the 
renormalization group coefficient f3(a s ), conformal symmetry is preserved in perturbative QCD. 
Even in the case of a non-vanishing (3 function the conformal representations can be changed in 
such a way that the evolution operator is diagonal. The evolution equation can then be solved 
trivially. The conformal moments at the input scale depend on the skewness parameter and can be 
expanded in an appropriate orthogonal polynomial basis where the expansion coefficients depend 
on the momentum transfer squared. In other words, GPDs and GDAs can be represented by 
a conformal partial wave expansion, where the expansion coefficients are characterized by form 
factors that are labelled by the conformal spin and by an appropriate second quantum number, 
e.g., the angular momentum. Most importantly, it has been demonstrated that the first few form 
factors are measurable on the lattice. Moreover, the crossing relation between GPDs and GDAs 
are reduced in this representation of the continuation of these form factors from the space- to the 
time-like region and reverse. 
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Group theoretical discussions based on discrete conformal spin have a long history in QCD. 
However, in practice, they seemed to be only useful for distribution amplitudes and GDAs, where 
the resulting series convergence [2U ESj- Combined with conformal symmetry predictions, the 
perturbative corrections for (virtual) two photon processes in the generalized Bjorken limit can be 
worked out, e.g., for the photon-to-pion transition form factor, to next-to-next-to-leading (NNLO) 
order accuracy [2E1 EZ| . 

For GPDs the conformal partial wave expansion, where the conformal spin is a non-negative 
integer, is represented by a series in terms of mathematical distributions, which only converges if 
it is convoluted with suitable test functions. An appropriate resummation of this series has been 
proposed in Ref. [2E] and the details are presented here. There are also several other suggestions 
in the literature to define out of this divergent series. One might insert the identity expanded in 
terms of polynomials, which is however only applicable for a certain kinematical region [2E1 ED] , 
and for a next-to-leading (NLO) analysis see Refs. [S3E21- One can also represent the identity by 
its Fourier transform which makes contact to the group theoretical representation with complex 
valued conformal spin [H3], later adopted for GPDs in Ref. [SI] and more recently in Ref. [3*H] . 
Also a resummation by an integral transformation has been suggested [SB] , which, however, under 
close scrutiny turned out to be unpracticable or at least rather complicated [37]. An attempt to 
approximately resum the conformal partial wave expansion within a Taylor expansion of conformal 
moments has been suggested in Ref. [HEj- Of course, all these proposals can be related to each 
other. However, because of the intricate mathematics involved a correct, complete and efficient 
resummation of conformal partial waves has not yet been worked out. 

In this paper we employ the Sommerfeld- Watson transformation to resum the conformal par- 
tial wave expansion of GPDs, finally yielding a Mellin-Barnes integral for GPDs. The resulting 
representation is similar to the one recently proposed in Ref. [32], however, not identical. The 
Sommerfeld- Watson transformation requires the analytic continuation of the conformal spin, which 
plays here the analogous role to the complex spin J in the inverse Mellin representation of parton 
densities. While the analytic continuation of Mellin moments for parton distribution functions 
is a rather simple task, it is a highly non-trivial one for the conformal moments of GPDs. This 
central problem is solved by us to an extent such that the framework can be applied as soon as a 
GPD ansatz is given. Although the final GPD representation as a Mellin-Barnes integral over the 
complex conformal spin seems to be rather complicated it has several advantages. The dependence 
on kinematic variables is separated in this representation, it allows a simple and stable numerical 
treatment of GPDs and their convolution with hard-scattering amplitudes, and can be used for 
the analytic approximation of scattering amplitudes. Moreover, the evolution equations to leading 
order (LO) accuracy are trivially solved and the conformal approach in Ref. |2oH2T| can be adopted 
for the study of higher order corrections in perturbative QCD. Especially, the NNLO corrections 
to deeply virtual Compton scattering (DVCS) are calculable in a rather economic manner [3*H] . 
Also the use of the crossing relation between GDAs and GPDs is possible in our representation 
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and should be very valuable for phenomenology. 

The paper is organized as following. Sect. is devoted to the conformal partial wave decom- 
position of GPDs for complex conformal spin. We start with a review on the anatomy of GPDs 
in Sect. 12. ll and discuss the extension of the GPD support [2E]- To the best of our knowledge this 
issue has not been presented in detail so far. Here it will guide us to find the correct treatment of 
partial waves with complex conformal spin. We consider then the crossing relation between GDAs 
and GPDs and derive in Sect. 12 .HI the new GPD and GDA representations in terms of Mellin- 
Barnes integrals. In Sect. 01 we consider evolution kernels and their convolution with GPDs in a 
scheme that preserves conformal symmetry. Moreover, we present the Mellin-Barnes integral for 
the scattering amplitude in DVCS and discuss its analytic approximation. In Sect. |I]we have a 
closer look to the analytic continuation procedure of conformal moments for a simple GPD toy 
ansatz. Then we address the issue of ansatze for conformal moments and explore the features of 
the resulting GPDs and GDAs. Furthermore, for vanishing longitudinal momentum fraction in the 
t channel we have a short look at valence quark GPDs. Motivated by lattice results [2B 122 122] , 
we introduce a parameterization for which the experimental constraints on GPDs indicate that 
leading Regge trajectories are present in conformal moments. Finally, we summarize and con- 
clude. In Appendix El integrals, which are used in the main text, and the rotation from ordinary 
Mellin moments to conformal ones are presented. Appendices [B] and contain the Mellin-Barnes 
integrals for gluonic GPDs and conformal evolution kernels, respectively. 



2 Features and parameterization of GPDs 
2.1 The anatomy of GPDs 

GPDs are defined as Fourier transform of light-ray operators, sandwiched between the initial and 
final hadronic states. There is a whole compendium of GPDs for each hadron. In addition the 
initial and final states can have different quantum numbers (transition GPDs), and one even can 
replace the hadrons by nuclei (nucleus GPDs). Once the initial and final states are specified, 
GPDs are classified with respect to the twist of the operators and the spin content of fields. At 
leading twist-two level three different types of quark and gluon GPDs can be defined (here the 
gauge link is omitted): 

9FV ) M f 7+ 1 

W A Ux, V ,A*,^) = /^e^+<P 2j 5' 3 |^(-icn)^7 + 7 i5 U;W|i'i^ 1 >, W 
W T ) 71 (ia +± ) 

g F A \(x jVj A\^)=2 J e iKxP+ (P 2 , S 2 \G a +ll (-Kn)l ie^+ X?2 + (kti)|Pi, S x ), (2) 

F I I Tinj'nfT \ 
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with P + = n ■ (Pi + P 2 ), V- = n* ■ V, n 2 = (n*) 2 — 0, n ■ n* — 1. In the first (vector) and second 
(axial- vector) entry the in- and outgoing partons have the same helicities, and the sum (vector) and 
difference (axial- vector) of left- and right-handed partons is taken, respectively. For the third entry, 
called transversity, a helicity flip appears. GPDs depend on the momentum fraction x, conjugated 
to the light-cone distance 2k, the longitudinal momentum fraction rj = (Pi — P 2 ) + / (Pi + P 2 ) + 
in the t-channel 2 , the momentum transfer A 2 = t = (P 2 — -Pi) 2 , and the renormalization scale 
fi 2 . The latter is induced by the renormalization prescription of the operators, which is part of 
the GPD definition. To deal with the polarization of the hadronic states, one might introduces a 
form factor decomposition j¥01 E] • For instance, for the nucleon GPD Dirac and Pauli-like form 
factors appear in the vector case [30] : 

l F V = U(P 2 , S 2 ) 1+ U(P U Si)H t ( X) rj, A 2 ) + U(P 2 , S 2 f-^^U(P U Si)E t ( X) r), A 2 ) , (3) 

where i = u,d, s, ■ ■ ■ ,G. To avoid confusion, let us note that for the process 7*(q , i) + p(Pi) — * 
7* (q 2 ) + p(P 2 ) two scaling variables exist. They are denoted as £ and t] and are defined by [2]: 

— q 2 A • q 1 

£ = "5 — ' ^ = ~~5 — with ^ = o(?l + ?2). (4) 
P ■ q P ■ q 2 

Both variable coincide, up to power suppressed corrections 0(A 2 /Q 2 ), when the outgoing photon 
is real, i.e., for DVCS one can simply replace t] by £. In this paper we will treat the general case. 
The definitions (HJ)-© imply the basic properties of GPDs: 

• In the forward limit A— >0 helicity non-flip GPDs reduce to parton densities El HI IS] , e.g., 

Qi(x^ 2 ) = hm Hi(x,rj, A 2 ,/i 2 ) (5) 
and the helicity flip GPDs Ei decouple, but 

limE^^AWC (6) 

• The /^-dependence is governed by linear evolution equations [HE], which can be derived 
from the renormalization group equation of the light-ray operators (121 EH] • 

• Hermiticity [2] together with time reversal invariance jlU] leads to a definite symmetry with 
respect to the skewness parameter rj, e.g., Hi(x,rf) = Hi(x, —rj). 

• The Mellin moments of GPDs are expectation values of local twist-two operators: 

J dx x n V v (x, V , A 2 , Q 2 ) = -^n" ■ ■ ■ (P 2 , S 2 \S ^ qlfl0 ■ ■ -iD, n V q \Pu Si) , (7) 



2 In the literature 77 is now denoted as £, which also is the Bjorken like scaling variable in hard inelastic exclusive 
processes. To be precise, we distinguish between both variables. The sign convention for £ is fixed, for rj it is 
changing. For quantities which are even under reflection, i.e., rj — > —77, the sign convention is irrelevant. Here we 
define 77 in such a way that it corresponds to the variable £, commonly used in the definition of GPDs, too. 
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where D^=D^ — is the covariant derivative, acting as indicated by the arrows, and the 
operator S symmetrizes all indices and subtracts the traces. Lorentz covariance enforces 
that this moments are polynomials in 77. 

Furthermore, GPDs are constrained in the region x > [77 ] by the positivity of the norm in the 
Hilbert space of states. The most general form of such positivity bounds [HJ EHl EH] , known so far, 
are given as an infinite set of constraints (13 EH]- Such constraints can be alternatively understood 
within the representation of GPDs as overlap of light-cone wave functions |49j . 

Let us consider the support of a GPD in more detail 3 . A generic quark GPD F(x, 77, A 2 ), e.g., in 
the vector case, is related to a double distribution (DD) D(y, z, A 2 ) by the integral transformation 

/! r l ~\y\ 
dy dzx p 5(x-y-r]z) £>(|/,z,A 2 ), p = {0, 1} . (8) 
-1 J-l+\y\ 

Here D(y, z, A 2 ) is an even function in z so that F(x, rj, A 2 ) has the proper symmetric behavior 
under the exchange rj — > —rj. Obviously, its Mellin moments, i.e., Jdxx n F(x, 77, A 2 ) are even 
polynomials in 77, since the support of D(y, z, A 2 ) is restricted. Depending on the form factors 
appearing in the decomposition of the GPDs (JTJ) and (J2J), see, for instance, Eq. 0, the order of the 
polynomial is n or n + 1 . To treat both cases in a convenient and generic manner, we have included 
in Eq. (JHJ) the factor x p with p = (p = 1) in the former (latter) case [H2]- This restores the 
correct order of the polynomials 4 . We can now fix 77 to be positive and decompose the integration 
with respect to y into y > and y < 0. This results into a decomposition of F(x, 77, A 2 ) in its 
quark q and anti-quark q part: 

F(x, 77, A 2 ) = q(x, 77, A 2 ) t q(~x, 77, A 2 ) . (9) 

Here both functions separately satisfy the polynomiality condition: 

q(x,r],A 2 )= dy dz x p 5(x - y - r] z) D(y,z,A 2 ), p = {0, 1} (10) 

JO J-l+y 

and analogous for the anti-quark GPD q(x, 77, A 2 ), where D(y, z, A 2 ) is replaced by D(—y, z, A 2 ) = 
±(-l)PD(-y,z,A 2 ) 

q(x,r ] ,A 2 )= [ dy [ " dz x p 6(x - y - V z) D(y, z, A 2 ) , p = {0, 1} . (11) 

J0 J-l+y 

3 The GPD support might be directly derived by means of a partonic Fock state decomposition and the so-called 
a-representation for Feynman diagrams, see for instance Ref. [5()j . Equivalently, a GPD can be expressed in terms 
of a DD pil51j. which has a simpler structure, and one might consider it as more convenient to derive the support 
of the former from that of the latter. 

4 Note that within p = an additive so-called D-term was proposed to generate r] n+1 terms - It is only 
non-zero in the restricted region \x\ < \rj\ and is contained in our representation with p = 1 as an additive term of 
.D(y, z, A 2 ) that is proportional to 5(y). Our parameterization offers the possibility that the rj n+1 terms arise from 
an uniform GPD. 
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X\PyS X2P1 



X\P\/ \x 2 Pi 



Figure 1: Partonic interpretation of GPDs in the central (left) and outer (right) region. 

Obviously, both quark and anti-quark GPDs have the same mathematical representation and so 
we will in the following mainly deal with the quark one. The results for anti-quark GPDs are 
easily obtained by replacements D — > D. 

In the central (exclusive or ER-BL) region —77 < x < 77, q(x, 77, A 2 ) might be interpreted 
as probability amplitude to have a meson like configuration inside the hadron, while the outer 
(inclusive or DGLAP) region 77 < x < 1 can be viewed as probability amplitude for emission and 
absorbtion of a quark with momentum fraction X\P\ = j^-Pi and x<iP\ = j^p-Pi, respectively, 
see Fig. ^ Remarkably, both regions have a dual interpretation, namely, as meson and parton 
exchange in the t and s channel, respectively. 

Lorentz invariance ties both dual regions, which can be read off from the representation 
that ensures polynomiality. Suppose 77 > 0, the z integration in the double distribution represen- 
tation (fTUI) can be trivially performed 5 and leads to the support 



The GPD representation (|12|) is manifestly invariant under the transformation 77 — > —rj, especially, 
the support —rj < x < 1 remains untouched. 

In the central region the GPD is given by uj (x, 77, A 2 ) from which the outer region, determined 
by the symmetrized function uj (x, 77, A 2 ) + uj (x, —rj, A 2 ), can be restored, see Fig. |2l Let us have 
a closer look at this continuation. In the central region the integration variable in the integral 
(|13|) takes the values < y < < 1. The restriction of the second argument in the double 
distribution \z\ = \(x — y)/r]\ < 1 — y is ensured by the values of the lower and upper limit, see 
solid line in Fig. Efa)- The integration path, starting at y = and z = x/77, lies inside of the 
DD support as it must be. At the cross-over point x = 77 it starts at the support edge y = and 

5 Since 77 > 0, we have set S(x — y — zrj) = l/r)5(x/ri — y/rj — z) rather than to indicate the modulus 1/|7?|. The 
sign convention of —77) and its transformation under reflection 77 — > —77 avoids an overall sign(?7) factor in Eq. 
H12JI . This allows us to treat u>{x,rf) as a holomorphic function in the complex 77 plane. 



q(x, rj, A 2 ) = 9 (-77 < x < 1) u (x, 77, A 2 ) + 9 (77 < x < 1) u (x, -77, A 2 ) . 



(12) 



The function u follows from the y integration in Eq. ()10j) 




(13) 



7 
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-0.5 



(a) 



Figure 2: Support of a DD (the area surrounded by thick lines) and integration path for the 
calculation of u(x,rj), see Eq. for the central region x\ < t] (solid) and the outer region 

r\ < Xi (dashed) are shown in Fig. (a). The continuation of the integration path over the support 
boundary is indicated as dotted line. In Fig. (b) the support of the resulting GPD (|12|) is depicted. 



f— 2 rather than zero. We can now define an 

2^ 



z = 1 and in the outer region the lower limit is y 

(ambiguous) continuation of the DD support for \z\ > 1 — y by any smooth function D(y, z, A 2 ) 
symmetric in z. 



D(y,(x-y)/rj,A 2 



D(y,(x-y)/r),A 2 )d[y 



x 



V 



-D(y, (x-y)/ V ,A 2 )9 



x 



V 



y ■ (14) 



1 — rjj \1 — rj 

This provides the (ambiguous) continuation of uo (x, 77, A 2 ) into the outer region, u (x, —r], A 2 ) 
is obtained by reflection symmetry 77 — > —rj and adding both contributions leads to the integral 
representation 

1 



to 



(x,r/,A 2 ) +u(x,-r],A 2 ) 



1 + 77 



77 x-r, 

1-7, 



dyx p D(y, (x - y)/r), A 



(15) 



in which the integration runs only over the original support of the DD. Hence, the ambiguity in 
the continuation of u (x, 77, A 2 ) and uo (x, —r], A 2 ) drops out in their sum. 

Let us suppose that D(y, z, A 2 ) can be viewed as a holomorphic function of y and z inside its 
support and has branch cuts outside of it. Then the integration path in Fig. (J2J (a) can cross or 
go along such branch cuts. To deal with a unique definition of uj (x, rj) for all values of x we might 
define its value within its integral representation IjlHj) for rj < x by the principal value prescription 



1 < ■ s 1 < 
: oj{x + le, rj) + —u![x 



ie, rj) for 77 < x . 



2 v 1 " 2 

For illustration we give a simple example for the Radyushkin ansatz 



q(x, 77) 



dy 



1 V dzS(x-y- V z) p^-U(\z\/(l-y)) 

l+y 



i-y 



(16) 



(17) 
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where the A 2 dependence is disregarded. The parton density is parameterized by a toy ansatz 
l(y) = ~~ vY / B(a + 1, (3 + 1), which can be considered as building block for realistic param- 
eterizations. A popular ansatz for the profile function is 



U(z) 



U(z\b) 



'1-2 



2\b 



5(6 + 1,1/2)' 



B(x,y) 



r(x)r( y ) 

T(x + y) 



where the parameter b controls the strength of the skewness effect. This GPD ansatz can be 
evaluated in an analytically form in terms of hypergeometric functions for non-negative integer 
value of the parameter b. We might choose here for simplicity b = 0, i.e., 11(2,0) = 1/2. In the 
central region this toy GPD reads 



q(x,T}) 



Y{2 + a + (3) fx + r) 



27 1 T(2 + a)T(l + (3) 



V 



l+a 



l + a, 1 -p 
2 + a 



x + rj 
1 + 7] 



for 



7] < X < 7] . (19) 



The extension of the DD support corresponds to the analytic continuation of uj(x, rj) into the outer 
region and so we find: 



q(x,r)) 



Y{2 + a + (3) 
2T]T(2 + a)T(l + (3) 



X + T] 

1 + V 

X — 7] 
1 — 7] 

We remark that the rescaled distribution 
X 1 



l+a 



l+a 



1 + a, 1 - (3 

2 + a 
l+a,l-,5 

2 + a 



X + 7] 
1+7] 
X — 7] 



(20) 



V 



for x > 7] . 



Y p ~ l q 



y ' y ' 



A 



dy 



1-2/ 



dzX*6(X-yY-z) D(y,z,A 2 ), p={0,l} 



(21) 



l+y 



has technically the same support as it appears in evolution kernels. Defining w by 



vj(X, y, A' 



l+X 

i+y 



dyXW(y,X-Yy,A 2 ), p={0,l} 



[22] 



we have 



Y^q[y,^,A 2 ) =sign(l + y)^ 



X 



Y 



Y -X 

i + y 



w (x, y a 2 ) + 



X 
Y 



-X 
-Y 



(23) 



The evolution kernels have the form (}2"3*|) in the collinear limit A 2 = 0. This enables us to adopt 
results for the representation of GPDs to the representation of GDAs (see next section) and of 
kernels. Needless to say, the GPD ()12|) follows from X = x/t] and Y = I/77, where uj(x,7], A 2 ), see 
Eq. (|13j). is related to w(X, Y, A 2 ) by the formula 



to 



(x,r / ,A 2 )=r / ^ f-,-,A 2 V 

\7] 7] J 



(24) 
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2.2 Generalized distribution amplitudes and crossing 

Another ingredient we need for our derivation of a Mellin-Barnes representation of GPDs is their 
extension to 1 < 77. This immediately leads us to further non-perturbative distributions, the so- 
called GDAs j2ini, which are related to GPDs by crossing The GDAs, denoted as $(z, (, W 2 ), 
are defined in analogy to the GPDs in Eqs. Q and ((2j), however, the initial state is replaced by 
the vacuum and the final one contains two hadrons. For instance, the crossing analog of q F for 
a spin-zero target reads 

= |^e^ 1 - 2 ^+(0|^(-«:n) 7+ ^M|P 1 ,P 2 ) (25) 

with P + = n- (Pi + P 2 ). Here < z < 1 and 1 — z are the momentum fractions of the quark 
and anti-quark, respectively, which produce the hadron pair with invariant mass squared W 2 . 
< C < 1 is the momentum fraction of one of the hadrons. 

In the following we give a rather generic discussion of the crossing relation to GPDs, which 
neglects details about form factor decomposition or quantum numbers. We consider also only a 
quark GPD and its analog. For anti-quarks the only additional aspect, that one has to imple- 
ment, is the sign convention in the decomposition (JHJ). The GDA reads in terms of the rescaled 
distribution (J2IJ) with X = 1 - 2z, Y = 1 - 2(, and A 2 = W 2 as 



$(zX,W 2 ) = yp-i g _ _ A' (26) 

\Y Y J X=l-2z,Y=l-2(,A 2 =W 2 

Consequently, from Eq. we read off its representation 

$(2, C, W 2 ) = 9 (z - C) w (1 - 2z, 1 - 2C, W 2 ) + 9 (C - z) w (2z - 1, 2( - 1, W 2 ) , (27) 

where the function w is given by the integral (|22|) . Having in mind that in Eq. f|26|) a rescaled 
GPD appears, it remains a trivial exercise to directly relate GPDs and GDAs 6 : 

<b(z,(,W 2 ), I fii'-'^ii.i')! , „. . 1 , „, . 1 „, 2 



.(l-^l-Wj l^ w(l ,,,A>)j' l-2,«-,l-2C«-,^«A^,(28) 

where the lower line is in fact the relation ((23$ . 

Let us shortly discuss this crossing relation. A GDA is obtained from the GPD analog by 

w (1 - 2s, 1 - 2C, W 2 ) = (1 - 2Cr 1 w H^ 2 ) • (29) 

Since < Q < 1, the second argument of uj covers the region |l/(2£ — 1)| > 1. Thus, the crossing 
relation requires to enter a kinematic region which is unphysical for GPDs, except for the point 



6 Note that the upper line in Eq. I)28|) is only valid for 1 < 77, i.e., ( < 1/2. For 1/2 < £ one must replace i] 1 ^ 10 
by sign(77)?7 1_p . It is more convenient to use the second line, valid for all values of 77, together with the explicit 
representations l|12|) and (|27|l . see also footnote [5J 
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rj — 1, i.e., C = 0. However, for a given functional form of a GPD, the relation (|29|) allows to fix 
its phenomenological parameters. So for instance, the same function appears after factorization 
in hard exclusive electroproduction of a photon or mesons on a hadron target as GPD and in the 
production of a hadron pair due to two photon fusion as GDA. The knowledge of the analytic 
form of the GPD in the central region is sufficient to perform the symmetry transformation (J2*§|) to 
obtain the corresponding GDA. Reversely, a GPD follows from a given GDA using the symmetry 
transformation (|24j). As expected from general reasons, the physical and unphysical regions are 
again connected by crossing. Moreover, as we realized above, u(x,r), A 2 ) is not uniquely defined 
in the outer region. As explained above this problem is artificial and does not affect the net 
contribution in this region. 

Let us stress ones more that the extension procedure of the support is unique, which was 
shown in connection with the support extension of evolution kernels E] . Here we adopt the 
same arguments. Suppose we know the function w (1 — 2z, 1 — 2(, A 2 ) in the region < £ < z < 1, 
which is equivalently to the knowledge of the GDA $(z, £, W 2 ), see Eq. (|27|). Next representing 
the GDA (|2T)j) in terms of the DD (j21jl . we realize that its convolution with any holomorphic test 
function t(z) yields an holomorphic function in £. For instance, one finds for p = 

£dzr{z)*{z, (, W 2 ) = \ J\y J *J dzr (^~ y{1 ~ 2C) ~ D(y,z,W 2 ). (30) 

Hence, also the Fourier transform of the GDA with respect to z is a holomorphic function in the 
conjugate variables A and £■ Consequently, we can employ analytic continuation. Then the inverse 
Fourier transform together with the crossing relation (|2*Hj) yields the result we desire, 

—e~ iXx ^AC 2 / dze iXil - 2z) <S>{z,(,W 2 ) , (31) 
■oo2vr L 

where ( = (rj — l)/2r/, W 2 = A 2 and AC denotes the analytic continuation of both variables A 
and £. (Actually, the variables stays real and analytic continuation is only used to extend their 
numerical values on the real axis, e.g., for A up to ± infinity.) 



2.3 Complex collinear conformal spin partial wave expansion 

We derive now a new representation for GPDs which has several advantages, already mentioned in 
the introduction. In fact we will deal with a partial wave decomposition of GPDs, where the partial 
waves are labelled by the complex conformal spin, the quantum number which characterize the 
multiplets (towers) of conformal operators. This is rather analogous to the partial wave expansion 
of scattering amplitudes with respect to the complex angular momentum, however, requires a 
more attentive consideration. Irrespectively, of whether the symmetry is preserved or not, one can 
introduce such an expansion. Certainly, conformal symmetry is broken in the non-perturbative 
QCD sector. Fortunately, up to calculable corrections proportional the non-vanishing (3 function, 
it holds true in the perturbative sector and thus has nevertheless still predictive power. 
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In the following we consider only quark GPDs, the gluon case can be treated analogously and 
the results are collected in Appendix [BJ To derive the partial wave decomposition of GPDs, the 
following steps will be performed: 

• First the conformal moments and partial wave expansion of GPDs are introduced for discrete 
conformal spin. Here GPDs are represented as a divergent series in terms of mathematical 
distributions. 

• This series will be summed in the unphysical region, where it can be viewed as an ordi- 
nary expansion in terms of orthogonal polynomials, by means of the Sommerfeld-Watson 
transformation, which requires the analytic continuation of the conformal spin. 

• Finally, we complete the Sommerfeld-Watson transformation and derive a representation of 
GPDs in terms of a Mellin-Barnes integral for the central region. The outer region follows 
from an suitable continuation in x that arises from the analytic structure of GPDs. 

2.3.1 Conformal partial wave expansion 

So-called conformal moments of quark GPDs are formed with respect to Gegenbauer polynomials 
r) n Cn^ 2 (x/rj) with index 3/2 and order n. These moments are given by the expectation value 
of local conformal operators, see below Eq. (|37|). Relevant group theoretical aspects can be 
found in Appendix [B] and in the review jHS] . These moments can be viewed as an appropriate 
generalization of the ordinary forward Mellin moments, used in the analysis of deep inelastic 
scattering. The Gegenbauer polynomials are orthogonal polynomials, possess a definite reflection 
symmetry Cn 2 (x) = (— l) n Cn 2 (— x), and are the only solution of the second order differential 
equation 

dx 

that is finite at the singular points x = ±1. These polynomials form a complete basis in the 
interval [—1, 1]. Here we rescale the polynomials and choose the normalization 

c n (x,v)=V n c n $ with <x) = T fJ^^^ Cll\x) (33) 

in such a way that in the forward case the ordinary Mellin moments appear: 

lim Cn(x, rj) — x n . (34) 

There are several possibilities to express c n (x, rj) in terms of hypergeometric functions, which might 
provide different prescriptions for the analytic continuation of the discrete variable n. Below we 
will use 



, v l - x 2 )Cf(x) = -{n + l)(n + 2)C 3 J 2 (x) (32) 



r(3/2)r(3 + j) p f-j,j + 3 



2 1 +JT(3/2 + j) ' z 1 V 2 



Ci&V)= almoin A WA o hzr) (35) 



rj — x 



2'/ 
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for complex valued j. Equivalently, we can express it in terms of (associated) Legendre functions 
of the first kind. 

The conformal moments of a GPD, separated into quark and anti-quark ones, are defined as 
m n (r],A 2 )= dxc n (x,r])q(x,r],A 2 ) , m n (r],A 2 ) = dx c n (x, r))q(x, rj, A 2 ) . (36) 
They are given by the expectation values of collinear conformal operators, e.g., in the vector case, 

m n ( Vl A 2 ) - (-l) n m n (?7, A 2 ) = ^|y^y^ (A, ^7+ C*<* fe±) V q \Pi,S 1 ) . (37) 

The rotation to the ordinary Mellin moments (J7J) and its inversion are given in Appendix [X] by 
Eqs. (|169|) and (|17(J|) . The operators are characterized by the conformal spin, which in our case is 
n + 2. The (conformal) moments can be either calculated on the lattice 7 , can be directly modelled, 
or evaluated from a given GPD ansatz. In the latter case they are naturally decomposed as, e.g., 
for quarks, 

m n (r], A 2 ) = fi n (r], A 2 ) + fi n {-f], A 2 ) , fi n (r],A 2 )= dxc n (x,r]) uj(x,r], A 2 ) . (38) 

J-ri 

This is a simple consequence of the symmetry relations c n (x,r]) = c n (x,—r]) together with the 
representation (fT^j) . The /i n (?7, A 2 ) are only defined in terms of u(x, r], A 2 ) and, thus, they can be 
quite general functions of rj. However, after symmetrization with respect to rj, see first formula in 
Eq. (}38|) . one obtains the polynomial m n (rj, A 2 ). 

Now we would like to invert the transformation (J36|) . As mentioned above the polynomials 
c n (x, rj), see Eqs. (J3l?j) and (jMj) . form only a complete basis in the central region [—77,77]. Let us 
denote by p n (x,r]) the polynomials that include the weight (1 — x 2 ) and an appropriate normal- 
ization 

Pn(x,v)=V* 1 Pn(jj) , Pn{x) = 0(1 - \ x \ ^(3/2)F(3 + n) t 1 " ^ Cl ' 2{ - X) ■ (39) 

The orthogonality relation for Gegenbauer polynomials reads in our notation 

1 

dxp n (x, rj)c m (x, 77) = {-l) n 5 nm . (40) 

1 

The minus sign in the argument of Gegenbauer polynomials in Eq. (}3T)j) is conventionally and 
induces the factor (— l) n in the orthogonality relation. This sign convention is appropriate to 
perform the steps which follow. Note also that the support restriction is explicitly contained in 

7 The separation of valence quarks and sea quarks is done by measuring even and odd moments. For instance, in 
the vector case: if the quark and anti-quark sea is equivalent, in even moments the complete sea drops out, while 
for odd moments valence and sea quarks are added. 
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the definition (|39|) and so the integration region in the integral (J4*U|) is restricted to the central 
region. We might now expand a GPD in terms of such polynomials (|39j) 

oo 

q(x, n, A 2 ) = ^(-l)> n (x, rj)m n (rj, A 2 ) . (41) 

n=0 

It is easily to see that the conformal moments (|36|) of a GPD are reproduced by this series (}4~Tj) . 
However, it is divergent as expansion in terms of polyonomials 8 . Especially, the restricted support 
property of each individual term does not imply that the GPD vanishes in the outer region. 
Rather one should understand this expansion as an ill-convergent sum of distributions (in the 
mathematical sense) that yields a result which is non-zero in the outer region. Indeed, p n (x,r]) 
can be considered as the nth derivative of a smeared 5 function: 

" = Jm0h, Sj u (1 - 1,2) " +li< " >( * - • (42) 

Taking now the forward limit \im v _+ p n (x , 77) = (x) / n\ , the series (f4*Tj) turns out to be the 
expansion of parton densities in terms of derivatives of 5(x). 

On the other hand, Eq. ()41|) might be a convergent series for 77 > 1 and by means of the 
crossing relation (J28|) we find for a GDA the partial wave decomposition: 

00 

C, W 2 ) = ]T(-1)>„(1 - 2z, l)M n (C, W 2 ) , (43) 

n=0 

where the polynomials M n ((, W 2 ) = (l — 2() n m n (l/(l — 2(), W 2 ) are of order n. If $(2, C, W 2 ) is a 
smooth function that vanishes at the end-points z = {0, 1}, the series converges and the moments 
M n (C, W 2 ) behave for n — »• 00 as ~ 2~ n n~ l l 2 ~ e with e > . Here the exponential suppression by 
2~ n is a consequence of our normalization, which has been adopted from the Mellin moments of 
GPDs respectively parton densities. 

The series (|41|) and (J43|) are the conformal partial wave expansions with respect to the con- 
formal spin n + 2. They have the advantage that in perturbative QCD the conformal spin is, to 
some extend, a good quantum number. So for instance the LO evolution kernels are diagonal with 
respect to Gegenbauer polynomials. Furthermore, in this expansion the x and A 2 (respectively 
z and W 2 ) dependence factorizes. In the case of GPDs the (x,r]) dependence is decomposed in 
an intrinsic x/rj and a remaining (77, A 2 ) dependence contained in the conformal moments. The 
(77, A 2 ) or ((, W 2 ) dependencies can be separated by an expansion of the moments with respect to 
an appropriate set of orthogonal polynomials. It has been proposed to expand M n ((, W 2 ) in terms 
of Legendre polynomials P/(cos(#)) [HE], t ne eigenfunctions of the rotation group SO (3) for spin- 
less states. Here cos(#) ~1 — 2^ = 1/77 and 9 is the scattering angle in the center-of-mass system. 
Certainly, it is appealing to have such an expansion with respect to angular momentum /, which 

8 This is also the case in the central region [ — 77, 77] , since the coefficients in front or the Gegenbauer polynomials 
are enhanced by the factor r\~^~ x , which divergences for |?/| < 1 at j — ► 00. 
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Y 



(a) (6) 

Figure 3: (a) The integration contour in Eq. (J4l)j) enclosing the real axis in the complex j-plane. 
Adding semicircles results in an integration path which is parallel to the imaginary axis, (b) The 
contour of the Schlafli integral ()46|) in the complex u-plane. 

makes contact to the SO (3) partial wave expansion. When spin is involved, the expansion is rather 
given in terms of Wigner functions, expressible in terms of associated Legendre polynomials. To 
be not specific on the spin content and in view of the definition in terms of conformal operators, it 
appears also natural to expand conformal moments with respect to Gegenbauer polynomials with 
index 3/2: 

n n 

m n (r), A 2 ) = Fnk(A 2 ) f] n c k {l/r,) , M n ((, W 2 ) = ^ F nk (W 2 ) c k (l - 2() , (44) 

fc=0 fe=0 

where all A 2 and W 2 dependence is absorbed into the form factors 9 F nk . Obviously, the crossing 
of GPDs and GDAs, i.e., the transfer from the space-like to the time-like region, concerns now 
only these form factors. The form factors, appearing in a specific partial wave expansion with 
respect to the angular momentum, are simply obtained by a rotation from the conformal ones, 
i.e., F nk (A 2 ). 

2.3.2 Sommerfeld- Watson transformation 

As already mentioned, the series ()41|) can not be directly used in practical calculation. Rather, 
it must be either resummed or the individual terms must be smeared by inserting the identity, 
expanded with respect to an appropriate basis. However, the latter method is only applicable in 
a restricted kinematical region and has been performed only approximately. So we chose to resum 
the conformal partial waves series for GPDs instead. We consider it first in the unphysical region 
7] > 1 and rewrite q(x, r], A 2 ) as a contour integral in the complex plane that includes the positive 

9 We note that in this expansion the scaling invariance and so the conformal symmetry is broken, since in 
general the form factors contain now a massive parameter for dimensional reason. This is not surprising, since in 
non-perturbative QCD this symmetry does not hold - otherwise there would exist only massless hadrons. 
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real axis, see Fig. El 

I r(oo) y 

q(x, 77, A 2 ) = — & dj , .s Pj{x,rj) mj(rj, A 2 ) . (45) 
2z y (0) sm(Trj) 

Here we included a factor l/sin(7rj), which has the residue ReSj =n 1/ sin (717) = (— l) n /7r for 
n = 0, 1, 2, • • -. Thus, if no other singularities are present inside the integration contour, the 
residue theorem leads to the conformal partial wave expansion (j41|) . The main difficulty is to find 
an appropriate analytic continuation 10 of both functions pj(x,rj) and rrij(r], A 2 ) with respect to 
the conformal spin n + 2. 

First let us define the analytic continuation of pj(x,rj). This can be done using its definition 
(j39jl in terms of hypergeometric functions. To include also the support restriction, we represent 
the analytic continuation of the Gegenbauer polynomials by the Schlafli integral 

, N r(5/2 + j) 1 , (u 2 -lV +1 
Pi(x,ri) = - — ) '—, JJ r & du) V— . 46 

pa i) r(i/2)r(2 + j)27T2 7 ( _ 1+e) (x + u V y+ 1 1 ; 

Here the integration contour is essentially the unit circle in the complex w-plane where the points 
— 1 and +1 are included, see Fig. |3] (b). The integrand has four branch points in the complex 
w-plane, namely at {— oo, — 1, — x/r], 1}. These points will be connected by a single branch cut 
that goes along the real axis from — oo to Max(— x/r], 1). It is easy to see that for non-negative 
integer j = n the Schlafli integral is equivalent to the definition (|39|). The integrand possesses 
now only a pole of order n + 1 at u = —x/t] (and at infinity). For |a;| < t] the pole is inside the 
integration contour and the residue theorem gives p n (x, rj). On the other hand for \x\ > r\ the pole 
is moved out of the contour and so the integral vanishes. 

Before we evaluate the integral ()46j) in general, let us consider the forward case r\ — 0. Here 
the integrand is essentially reduced to the function (u 2 — l) jf+1 and possesses for non- integer j a 
discontinuity on the real axis in the interval — 1 < u < 1. We might now deform the contour so 
that the real axis is pinched. For |3?eM| < 1 we pick up a phase factor e ±J7r ^ +1 - ) for 3m« ^ and 
so we can write 



pJx, V = 0) = x^ 1 F(5 ( 2 / +J) r — (e^ J+1 ) - e" iir fdu (1 - u 2 y +1 . 

Fjy ' ' ' r(i/2)r(2 + j)27rz v J J-i 

The remaining integral represents just T(l/2)r(2 + j)/T(5/2 + j), which results in 

p^ V = 0) = Sin ^ [j + 1]) x-i-K (48) 

7T 

This is, up to the conventional factor sin(7r[l + j])/7r, nothing else but the integral kernel of the 
inverse Mellin transform, widely used in deep inelastic scattering. 

10 We remind that the analytic continuation of a function that depends on a discrete variable is not unique. For 
instance, a term proportional to sin(7r/) could be added, which drops out for j = n = 0, 1, 2, • • •. 
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For f) > and non integer values j the integrand in Eq. (ffify has no discontinuity for x < —77. 
For x > —J] we will pick up phase factors by surrounding the branch points at —x/r) and x — 1. 
They appear in the following interval above or below the real axis: 

-x/r), 1] 1 for j -77 < x < 77 and 77 > 
[-1,1] J I 77 <x 

Consequently, at the endpoint x = —rj and for x < —77 the integral vanishes 

< —77, 77) = . (50) 

For both the central and the non-vanishing outer region the complex integral can be evaluated by 
deforming the contour as before so that the real axis is pinched. Taking the discontinuity yield 
the following integrals along the real axis 

, . r(5/2 + j) sin(7rh' + l]) f 1 , (1 - u 2 V +1 



r(i/2)r(2 + J ) vr J^ xh (x + ur,y+i' 



r(5/2 + J ) sin(^ + l]) r du ^~^ (51) 



r(i/2)r(2 + j) 7r y_! (x + ur / y+ l 

At the cross over point x = rj both integrals have the same value and represent a Beta function. 
So pj(x = 77, 77) is smooth in the vicinity of this point and takes the value 

p (x - „ „) - 2l+J sin^b + l]) 

^(x-r;,^ - ^ + i r(3/2)r(3 + j) w • ^ 

The integrals in Eq. (|51|) define hypergeometric functions. The analytic continuation of the 
mathematical distributions p n (x,rj) with respect to the conformal spin is expressed by them as 
following: 

Pj (x,ri) = 9( v - IxDrj-^Vj (^) + 9{x - rfirj-^Qj (^) (53) 

where 



77/ \r) 



Here, a few comments are in order. First, for j = n = 0, 1, 2, ■ ■ • only the central region 
contributes and the relation 



2F 1 



2 



A 2 f-j-lj + 2 

1 — x \ 2 



x 



(56) 



establishes the definition (|3T)|) of p n in terms of Gegenbauer polynomials, cf. Eqs. (f3l?J) - (f33j) . Ob- 
viously, in the central region the analytic continuation is based on the definition of hypergeometric 



17 




Figure 4: The real (left) and imaginary (right) part of the conformal partial wave Pj(x,rj = 1) 
(solid), its first (dashed), and second (dotted) derivative for j = — 1/4 + i/2. 

functions. In the outer region, however, the result might be surprising. To clarify its meaning, we 
decompose the integral (Jo"Tj) for rj < x as J^du ■ ■ ■ = f^ x i du ■ ■ ■ — J_\du ■ ■ ■ and realize that it 
can be expressed by the function Vj. Thus, we can write Eq. (|53|) as 



Pj(x,rj) = 6(-r] < x) r]~' J ~ x V 3 (^j + 9(r) < x) cos(tt[j + 1]) rf^Vj 



(57) 



Here it is understood that the principal value is taken at the branch cut, which starts at x — rj, 
i.e., we insert [Pj(x + ie) + Vj(x — ie)]/2 for x > 1. The cos(7r[j + 1]) term in the second expression 
on the r.h.s. arise from the continuation of rj to —rj, again by taking the principal value. Hence, 
one realizes that this result precisely fits the structure of the representation (fT2j) for q(x, rj, A 2 ) in 
terms of the functions u>. We remark that the identity 

Qj(x) = ~Vj(x + ie) + ^Pj(x -ie) + cos(Tr\j + l])Vj(-x) for x>l, (58) 

we derived here, is a known relation between associated Legendre functions of the first and second 
kind. 

The function pj(x,rj) is continuous at the cross-over point x = rj, however, the imaginary part 
of the first derivative has a jump, while the real part is still continuous, see Fig. 0] It satisfies the 
second order differential equation for conformal partial waves with complex valued conformal spin 
j + 2: 

d 2 

(1 - x 2 )^ Pj (x,rj = 1) = -(j + l)(j + 2)p 3 {x,r, = 1) . (59) 

The expressions for the conformal partial waves in Eqs. (|53|). (|54j) . and (J55|) . turn out to be 
in agreement with a representation, given recently in Ref. (35], in terms of associated Legendre 
functions for the gluonic GPDs. In fact, up to some normalization factors and shift in both the 
conformal spin and the index, the same functions appear in the central and outer region. For the 
former one we, however, prefer to work with complex conformal spin, too. This incorporates the 
underlying duality between central and outer region in a manifest manner and yields a uniform 
representation of scattering amplitudes, see below Sect. 13.21 This is essential for the perturbative 
QCD analysis at higher orders. 
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The analytic continuation of the polynomials m n (rj, A 2 ) is denoted as rrij(r], A 2 ). These func- 
tions will be also analytic in 77, however, might have branch points at rj — 0, r\ — ±1, and rj = 00. 
It would be desirable to have an integral representation that makes this property transparent and 
might allow the continuation from 77 > 1 to 77 < 1 or even to negative values. Moreover, we 
will also require that the moments 771^(77, A 2 ) are bounded at large j. It turned out that with 
these requirements the analytic or numerical calculation of moments from a given GPD is a rather 
intricate task. We have to admit that this mathematical problem is also not solved here for 
any conceivable GPD. In the following we give, however, some recipes to evaluate the conformal 
moments for complex conformal spin in the region |?7| < 1. 

To derive an appropriate integral representation, we decompose here, in contrast to Eq. (jEEj) , 
the conformal moments into contributions that arise from the outer and the central region: 

m n ( V , A 2 ) = ^ n (v, A 2 ) + (V, A 2 ) • (60) 

Again polynomiality is only manifest for the sum but not for the separate terms on the r.h.s. The 
analytic continuation of /i° ut (77, A 2 ) is defined in terms of hypergeometric functions (1851) 

/^ Ut (77,A 2 )= / dx Cj (x,T]) [^(£,77, A 2 ) +^(2,-/7, A 2 )] . (61) 
Jr) 

From the asymptotics of the hypergeometric functions 

2 1 for 3 -> 00 , |arg(j)| < tt/2, (62) 

we can estimate the large j-behavior of conformal moments 



ji^fo, A 2 ) ~ f\x ( X + v/ ;f ^ M*, V , A 2 ) + U ( X , -v, A 2 )] 



(63) 



For < 7] < 1 the integrand and thus the conformal moments are exponentially suppressed for 
j — ► oo with |arg(j)| < ir/2. In the limit 77 = we arrive at the parton densities and get the power 
like suppression factor j~ p . 

The contribution from the central region reads for integer values n 

f i c : n (r ] ,A 2 )= [ dxc^x^u^tf). (64) 

The analytic continuation of the conformal spin, as done for the central region above, would 
yield terms proportional to sin(7rj) that exponentially grow at large Smj. The presence of such 
terms can also be read off from the integral ()16H) . in the appendix, together with our final GPD 
representation (fTTj) . given below. In the case that this integral can be evaluated in an analytic 
form for integer n these terms drop out and the analytic continuation can be performed by means 
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Figure 5: The integration contour in Eq. (|65j). 



of the substitution n — > j. We remark that the explicit knowledge of fi° nt (r], A 2 ) allows to restore 
the contribution from the central region by means of the polynomiality condition. Obviously, 
within such a procedure one might miss some terms that satisfy the polynomiality condition, but 
contribute only in the central region. Nevertheless, this procedure offers in principle a way to 
restore the polynomiality of a GPD that is only known in the outer region, e.g., from an ansatz 
within the overlap representation |49j . 

Another possibility for the analytic continuation of the conformal spin, appropriate for numer- 
ical calculations, arises from the definition of the conformal moments in the central region by the 
following contour integral (see FigEJ): 



n+l rl-e 

/O^, A 2 ) = "tt - r dz d n( z ) «>fa> V, A 2 ) . (65) 

l+e 



2i 

Here the function d n (z) is defined in the complex plane by the integral 



d.W-i-' / V-f )«■_('■'> . ( 66) 



7T (1 — Z 2 ) J_ x X 

It has a branch cut on the real axis in the interval [—1,1] and in addition single poles at x — 
{ — 1, 1}. Thus, we excluded in the integration contour of the integral (JHH)) the points x = { — 1, 1}. 
Obviously, inserting Eq. (J55j) into Eq. (J55j) yield by means of the residue theorem and taking the 
limit e — > to the original definition of the conformal moments (|64j). The function d n (z) can be 
expressed through Legendre functions of the second kind and an appropriate continuation of the 
conformal spin n is given in terms of a hypergeometric function 

_ 2-^- 2 r(l+j)r(3 + j) a-'- 1 f(j + l)/2, (j + 2)/2 J_\ 
j[X) ~ r(3/2 + j)r(5/2 + j) a: 2 -l *\ 5/2 + j x 2 J ' [ } 

The function dj(x) behaves at large x as x _J ~ 3 and has a pole at x — ±1. At large j it is bound 
in the whole region 1 < x < oo by 

-i-j 



d j {x,7i)~2-S?j—. (6? 
x z — 1 



Let us suppose that the function lj(xt), rf) is a regular function that has only a branch cut along the 
real axis [— oo, —1] and has a power like behavior at infinity. We can expand the upper and lower 
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part of the integration contour in Fig. El to semicircles encompassing the complete half-planes such 
that the real axis is pinched in the intervals [—00, —1] and [1, 00]. For a certain value of 5ie j, the 
contributions at infinity are negligible and so we pick up the discontinuity along the negative real 
axis and the pole at x = 1. This result serves now for the definition of the analytic continuation 
in j: 

Hf n (r],A 2 \a) = ar] J+1 / dxdj(-x) — [cu(xr] - ie,r), A 2 ) + u(xr] + ie,r), A 2 )] 
r,v+T(l/2)r(l+i) , _ 



+ UJ V(3/2 + ,-) "("■"■ A) - (69) 

Our choice of Jte j guarantees the convergence of the integral and, moreover, we assumed here that 
oj(x = —r],r)) vanishes. The choice of Jtej can be relaxed by introducing appropriate subtraction 
terms that improve the power behavior of Lu(r)x, rj) at x — 00. On the r.h.s. of Eq. (J69|) the last 
term arises from the pole at x = 1 and the first term contains a phase factor a = e ±MrJ ', that 
depends on the continuation procedure from positive to negative real valued x. To get rid of this 
factor we define for even and odd n the analytic continuation as 

ACfoA")-/ ^ ( \^ k = 1) J ton=<H*. (70) 
Pn w ' ; 1 /xf n (77,A 2 |a= -1) J I odd v ; 

Note that the appearance of the phase factor is strongly connected to the analytic properties of 
GPDs and is only absent if the GPD has no branch point at x = —77. 

Finally, we have for the analytic continuation of the conformal moments 

mf A 2 ) 1 f ^fa, A> = 1) \ reven 

m n (77,A)^<^ aox V=^. (r?,A) + 4 A9I _ S for n= ,, . (71) 



f d (r/,A 2 )J n w ' 7 [ ^ ccn (r/,A 2 |a = -1) J I odd 



It is obvious that m™ and m° dd lead only for even and odd integer values of j = n to polynomials, 
which are even in 77. The difference of even and odd moments arises only from the central region 
and is expressed by a function in 77 which does not degenerate into a polynomial for integer values 
of j = n. 

2.3.3 GPDs represented as Mellin-Barnes integral 

Now we like to change the integration contour in Eq. (|45jl . in such a way that it becomes parallel 
to the imaginary axis in the j-plane, i.e., extends from c — ioo to c + ioo. Here the constant c < 
is chosen such that all singularities contained in the conformal moments mj(r], A 2 ) are on the left 
hand side of the integration path. Let us assume for the moment that rrij(r), A 2 ) degenerates for 
all non-negative integer values of j to polynomials. As displayed in Fig. El (a), we add two quarters 
of circles in the first and fourth quadrant so that the integration contour includes the imaginary 
axis and is closed by the arc that includes the points c + zoo, 00, and c — ioo. 
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It remains to show that within our definition of conformal moments the contribution from the 
arc vanishes. Suppose that rj > 1 and that the variable x is rescaled by 77, i.e., x = rjX. Hence, 
we must study the behavior of pj(r)X,rj) for j — »• 00 in the interval \X\ < 1. The asymptotic 
expansion of the partial waves for large j with |arg(j)| < n/2 can be read off from the behavior 



of hypergeometric functions jo6j and is 

Pj(Xr], 7]) ~ [e j ™c™ s H-X±ie) ± • e -(j+3) arccosh(-X±i e )) for |X| < 1 . (72) 

For — 1 < X < 1 the function arccosh(— X±ie) has a monotonously increasing or decreasing imag- 
inary part that lies in the interval [0, tt] and [— 7r, 0] for the +ie and — ie prescription, respectively. 
Thus, within both prescriptions 

1 

Pj(v x ,v) m j(v^ 2 ) (73) 



sin(Trj) 



exponentially vanishes for —1 < X < 1 on the arc, specified above, as long as rrij(rj, A 2 ) behaves 
for rj > 1 as 

7^(77, A 2 ) ~ f^Y for (74) 



.2 

This remains true, if we replace pj(x,rj) by the symmetric and antisymmetric combinations 

pf cn (x, 77) = - [pj(x, r/) + Pj(-x, r/)] , p° dd (x, 77) = - \pj(x, rj) - Pj(-x, rj)} (75) 

and rn,j(rj, A 2 ) by the corresponding even and odd moments, respectively. Since p n (x,r)) = 
(—l) n p n (—x,r]) has definite symmetry, the residues of 

-^-pf cn (x,r/)m J cvcn (r/,A 2 ) and -±-pf d ( x , rj)mf d (rj, A 2 ) (76) 
sm(7rj) J J sm(7rj) J J 

only contribute for even and odd values of j = n, respectively, and, thus, the polynomiality is 
implemented in a manifest manner. 

Unfortunately, we did not give in the previous section a representation for mj(rj, A 2 ) that allows 
the analytic continuation in rj to the region rj > 1 and simultaneously satisfies the requirement 
that rrij(rj, A 2 ) fulfills the bound ((73}. This can be read off from Eq. ()63|). where obviously, we will 
pick up an additional phase. It seems that such a representation can not easily be found rather 
we encounter here similar problems as for the continuation to the region rj < 0. 

This difficulty can be avoided once we replace in m,j(rj, A 2 ) the variable rj by rj* with 77* < 1. 
It is sufficient to consider the cross-over point, since for all other x values inside the central region 
we will have an additional exponential suppression due to the phases of pj (77, 77), see Eq. (172}) . For 
X = 1 we find from Eq. (|52jl that Pj(rj, rj) behaves on the cross-over point as ~ (2/r/)^ sin(7rj) / \fj. 
To get rid of the exponential growth, induced by the real part of j, we might choose 77 > 2 and can 
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arrange in this way even an exponential suppression. If now the conformal moments mj(r)*, A 2 ) 
for given 77* do not grow faster than (rj/2y /y/j for j — > 00, the integral on the infinite arc does 
not contribute for x — 1, too. Thus, after appropriate scaling of x with 77, analytic continuation to 
the region < 77 < 1, and setting 77* = 77 we find for all values \x\ < 77 the following Mellin-Barnes 
integral representation for GPDs: 

1 r-c+ioo 

q(x, r], A 2 ) = - / dj —— Pj (x, 7/) 77^(77, A 2 ) . (77) 

2 Jc-ioo sm(Trj) 

Finally, we extend this integral into the outer region. For pj(x,r]) this is done by means of the 
definition (jo^j) . 

The Mellin-Barnes integral ([??]) can be used when the analytic continuation of even and odd 
moments leads to the same function 771,(77, A 2 ). If this is not the case, we separately introduce the 
Mellin-Barnes integral for the even and odd part of 

q(x, 77, A 2 ) = g cvcn (x, 77, A 2 ) + q odd (x, 77, A 2 ) (78) 

Since the representation ()77)) remains valid if we substitute pj(x,rj) by the symmetrized partial 
waves py cn (x, 77) and p° dd (x, 77), we have the same form of the Mellin-Barnes integral as before: 

• rc+ioo 1 

g even/odd (X) ^ A 2 ) = M dj -^L-p^^x, t/) mf ^ odd ( v , A 2 ) . (79) 

^ Jc-ioo sm l 7r JJ 

Bearing in mind that pj(x,r]) is set to zero for x < —rj, the extension of pj ven (x, 77) and p° dd (x, 77), 
defined in Eq. (fTH|) in terms of pj(x,r]) and pj(—x,r]), is consistently done by the procedure (jSHJ)- 
The original GPD (J78j) is restored by adding its even and odd parts: 

• pc+ioo 1 

q(x, t/, A 2 ) = - / dj —— [pf en (x, t/) mf -(77, A 2 ) + p° dd (x, 77) m° dd (7/, A 2 )] . (80) 
2 Jc-ioo sm(vrj) 

The extension into the outer region is not based on analytic continuation and the reader might 
wonder whether this prescription (p)3|) is indeed correct. Besides the arguments we gave in Sect. 
12.3.21 we provide next some further support. Let us first verify that the GPD representation (JHUJ) 
has the correct support property, i.e., it contributes only for —77 < x. As mentioned above, the 
difference of odd and even moments arise only from the central region and has thus the functional 
form (j)/2y +l fj(r), A 2 ), where fj(r],A 2 ) has a power-like behavior at j — > 00, see Eq. f!69|) . We 
expect that such a term does not contribute to the outer region. Indeed this can be read off from 
the Mellin-Barnes integral, which takes the form 

■ rc+ioo O-j-l 

see definitions (|53|) and For x > rj the term 2~i~ 1 Qj(x/rj)/ sin(7rj) contains an exponential 

damping factor (rj/xY for j — >• 00 with | arg(j)| < 7r/2, while the remaining factor in the integrand 
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is bounded. So we can close the integration contour by an infinite arc that includes the first and 
forth quadrant. Since the whole integrand is analytic in these two quadrants, Cauchy theorem gives 
zero. Note that Qj{x/rf)/ sin (717) contains no singularities on the real positive axis. Consequently, 
even and odd contributions are the same in the outer region and will add for x > rj, while they 
cancel for x < —rj. 

Alternatively, we can represent the GPD as 

1 rc+ioo ^ 

q(x,r],A 2 ) = - / dj—- — [6(\x\ < rj)pA-x,rj)mf{r), A 2 ) + pAx,r))mf(r), A 2 )] , (81) 

2 Jc-ioc sm(TTj) 

where = (m J evcn + m° dd )/2 and = (m™-m° dd )/2. Here the polynomiality is not manifest 
in the conformal moments, rather it is separately restored for mf n + m^ n and rn^ n+1 — rn^ n+1 . 

Let us next show that the Mellin-Barnes integral vanishes for \x\ > 1. Here the expression 
7]~^~ l Qj{x /rj) I sin(7rj) in the integrand vanishes on the infinite arc, due to the damping factor 
(l/x)- 7 , see Eqs. f!53|) and (|55|). As in the previous paragraph, we can close the integration contour 
and find that the integral gives zero. If we would have allowed an additional term proportional 
to Vj(—x,rj) in the extension of the support for x > rj, it would in general contribute in the 
unphysical region x > 1, too. 

Finally, the correctness of the GPD representations as Mellin-Barnes integral is deduced from 
the fact that for non-negative integers the conformal moments are reproduced. Employing the 
symmetry with respect to x — > —x, we find from the representation (|80[). e.g., for even moments, 

) /-oo ^ r-c+ico -y 

dxc n (x,r])q(x,r],A 2 )= dxc n (x,r])- dj , .. Pj(x, rj) mf cn (r], A 2 ) (82) 

J-n 2 Jc-ioo sm(Trj) 3 

for n — 0, 2, ■ • •. From the moments separately calculated in the outer and central region, given 
in Eqs. (jl59|) and (jl61|) . respectively, it follows 11 

f°° dx c n (x, rj)q(x, rj, A 2 ) = lira J- dj ( ' J^iML + JML) mJ^A 2 ) , (83) 

for n = 0, 2, • • •. The integration path can be parameterized by j = n + iX and making use of 
the identity 1/(A — ie) — 1/(A + ie) = 2iir5(X), the integral yields finally the conformal moment 
m^ ven (?7, A 2 ) for n = 0, 2, ■ • -. Analogously, the conformal moments m° dd (r/, A 2 ) arise for odd 
values of n. 

We complete this section with the Mellin-Barnes representation for GDAs, which follows now 
by crossing from Eq. (|77jl with the constrain z > (, which arises from —l/r)<x/r)<l. Employing 

n The x-integral over the outer region only exist for 3?ej > n. Since here the integrand has no singularities on 
the positive real axis, we can first shift the Mellin-Barnes integration path to the imaginary axis to the right. Then 
with 5Re j = n + e the x- integration is performed. On the other hand the x- integration over the central region must 
be performed within the original contour. It produces a sin(7rj) term, which removes the poles on the real axis and 
so we can shift the Mellin-Barnes integration path along the positive real axis, too. Both remaining integrals can 
then be combined in the limit that is indicated in Eq. (|83|l . 
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the symmetry transformation z — > 1 — z and £ — > 1 — £ we get the remaining contribution, see 
Eq. flZZD: 



i r c+ioc i r f z — > 1 — z 1 

* (2 ' c - ^ = 2 jL * 5wi r " ° B(1 " 2z ' imc w2) + 1 c - 1 - c I 



(84) 



where Mj(£, iy 2 ) is the analytic continuation of the conformal GPD moments, obtained by cross- 
ing, see below Eq. ()43|). If Mj(£,W 2 ) is bounded along the integration path for < £ < 1, the 
integral of the first term in the square brackets in this formula exist for all values of < z < 1. 
So we can drop both the restriction 9(z — () and the second term in the square brackets, obtained 
by symmetry. This latter term follows now from analytic continuation and so we arrive at the 
representation: 



§{zX,W 2 )= l - I dj——p J (l-2z,l)M j ((,W 2 ). (85) 

^ J c-ioo SltiyKJ) 

Suppose that Mj((,W 2 ) vanishes for j — > oo with |argj| < tt/2, it is straightforward to see that 
the conformal moments for non-negative integer conformal spin are reproduced. Employing Eq. 
f| 1 6 lj) with 77 = 1 and k = n = {0,1,2,- - •}, we can then close the contour so that it now encircles 
the positive real axis. The residue theorem yields then M n (£, W 2 ). 
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3 Evolution kernels and coefficient functions in a manifest 
conformal scheme 

To LO accuracy the evolution kernels and coefficient functions for exclusive processes respect 
conformal symmetry, which is the symmetry of the QCD Lagrangian at the classical level for 
massless quarks. Conformal symmetry find, in fact, many practical applications in QCD. It allows 
for instance to solve the mixing problem of light-ray operators, caused by renormalization, and 
by means of the conformal operator product expansion it predicts the hard-scattering amplitude. 
The breaking of conformal symmetry beyond LO has been studied with perturbative methods in 
great detail |HZl EH1 EH EH1 EM EI] • The main lesson from these studies is that the only physical 
contributions violating conformal symmetry are generated by the non-zero f3 function. All other 
terms violating conformal symmetry (in the perturbative sector) are mainly artifacts introduced by 
the standard renormalization/factorization prescription for the light-ray operators, which is based 
on some version of minimal subtraction within the dimensional regularization scheme. Within 
the standard scheme the conformal symmetry is separately broken in the evolution kernels and 
the hard-scattering amplitude for two-photon processes. However, it can be restored by a finite 
renormalization providing a scheme in which conformal symmetry is manifest, except for terms 
proportional to (3 that are induced by the trace anomaly of the energy momentum tensor. Thus, 
we will restrict us in the following two subsections not to LO accuracy. Rather, the results can be 
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(at least partly) applied to all orders of perturbation theory. Of course, terms proportional to (3 
need special consideration. A first discussion of this issue can be found in Ref. [27], see also Ref. 



3.1 Convolution of GPDs with conformal kernels 

In this section we consider the convolution 



p-K(-X)q(y,V,^ 2 ) 



of a GPD with a generic kernel K(x, y) that respects conformal symmetry. Here the integration 
region is defined by the combined restrictions for the GPDs and the kernels. The requirement 
of conformal symmetry means in fact that the kernel has for \x\, \y\ < 1 the following spectral 
representation 

oo 

K{x,y)=^2{-l) n p n (x)k n c n {y) for |z|,|y|<l, (87) 

n=0 

where k n are the eigenvalues of K(x,y) and the polynomials p n (x) and c n (y) are defined in Eqs. 
(j2HJ) and (|39|). respectively. Suppose that even and odd eigenvalues have the same sign and provide 
after analytic continuation with respect to the conformal spin the same holomorphic function kj 
so that it satisfies a bound for all values of j with arg(j) < ir/2. In fact the eigenvalues k n coincide 
with the Mellin moments of the DGLAP kernel, given by rational functions and harmonic sums, 
and so we know their analytic continuation. We remark that the support in the (x, y)-plane 
is defined in Eq. f|193|) and the extension to the full region is unique, see Appendix The 
convolution of a GPD confined to [—77, 1] with this kernel leads to a function that has also the 
support — 77 < x < 1. 

To proceed in the simplest possible manner let us extend the representation (|87|) to the whole 
region by the series 

-K (-, A = V(-l)> n (x, V )k n c n (y, 77) . (88) 

Again, here we understand that p n (x, 77) are mathematical distributions with the restricted support 
—77 < x < 7] for integer values of n, while c n (y, if) are the polynomials ()33|) . which can be extended 
to the whole y region. Taking the Mellin-Barnes integral (|77|) and the spectral representation ()88j1 , 
the convolution (J86|) leads to a divergent series 

00 

K <g> q(x, rj) = ^(-l)> n (x, r])k n m n (r], A 2 ) . (89) 

n=0 

Here we employed the integrals (jl59j) and (jlfil)) for discrete conformal spin, as has been explained 
at the end of Sect. 12.3.31 Since the analytic continuation of k n does not spoil our assumptions for 
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the Sommerfeld- Watson transform, we can proceed as in Sect. 12.3.31 As expected, for a conformal 
kernel the convolution 

1 rc+ioo y 

K ® q(x, 77) = - / dj ——pjix, 77) kj mj(rj, A 2 ) (90) 

2 Jc-ioo sm(7rj) 

is in the Mellin momentum space given by a multiplication of the conformal GPD moments with 
the corresponding eigenvalues. 

In those cases in which we must distinguish between even and odd eigenvalues of the kernel, 
we write K as 

~ K (-> ") = D- 1 )" [pT n (x,v)kT n Cn(y,v) + P° n dd (x,r ] )k° dd c n (y,r ] )] , (91) 

where p^ cn (x) and p^ d (x) are defined as in Eq. (|7H|) . and employ the GPD representation (JHDJ). 
Since even and odd Gegenbauer polynomials and partial waves have definite symmetry under 
reflection, they can not mix. Hence, the convolution leads to 

■ rc+ioo 1 

K <g> q(x, rj) = - dj — — — [pf cn (x, rj) fcf en mf cn (r/, A 2 ) + {even -> odd}] . (92) 

If the support of the GPD was in the interval [—77, 1] and the eigenvalues are different in the 
even and odd sector, the convolution certainly gives us a function that lives now in the whole 
region [—1,1]. The support extension is caused by the mixing with an anti-quark GPD. With 
^even _ fcY, _|_ a an( j ^odd _ ^e _ fcA ^ we can d ecom p 0se the convolution as 

K (g) q(x, rf) = 5q(x, rf) + 5q(—x, rf) . (93) 

Based on the representation (|81|) . we interpret the terms on the r.h.s. as contributions to a quark 
GPD, with -77 < x < 1, 

• rc+ioo j^E 

5g(x, 77) = - / dj . / [0(|z| < rf)pA-x, rj)mf{rj, A 2 ) + p,(ar, r})mf(rj, A 2 )] (94) 
2 Jc-ioo sm(vrj) 

and an anti-quark GPD, with — 1 < — x < 77, 

2 rc+ioo fcA 

6q(-x,r)) = - J dj gin / < v)Pj(x,r])mf(r],A 2 ) + ^(-2, r))mf(r), A 2 )] . (95) 

For the convolution of a conformal kernel with a GDA we would obtain the analogous results. 
The convolution is conventionally defined as 

r 1 

K®$(z, C, W 2 ) = / dyK(l-2z, 1 - 2y)$(y,(,W 2 ) . (96) 

Note that the change of variable (1 — 2y) — > y, induces a factor 2 in comparison to the defini- 
tion ()86|) . We can now represent the GDA by the convergent series ()43|) . and the convolution 
immediately leads to 

00 , 
K&Q(z, C, W 2 ) = ^(-l)> n (l - 2z, l)yM n (C, V^ 2 ) . (97) 

n=0 
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Within the Mellin-Barnes representation the solution of the evolution equation is a trivial task 
to LO accuracy, where the conformal symmetry is manifest in any scheme. The restoration of this 
symmetry is well understood to NLO and even the terms proportional to (3 can be diagonalized 
with respect to conformal partial waves. Relying on this symmetry and borrowing the eigenvalues 
of the evolution kernels from the Mellin moments of the DGLAP kernel from Ref. |d^ IM] one can 
even proceed to NNLO. 

The simplest example is given by the flavor non-singlet sector and LO accuracy. Here the 
evolution equation reads 

^q(x, V , A 2 , „ 2 ) = -^ 7 (°) ® q( x , v , A 2 , , (98) 



where the evolution kernel is 

7 (°) (*,?/) 



r^l \ 1 + X f -i 2 \ r^l A - X 

®{x,y)-— 1 + +Q(-x,-y) 



(99) 



y \ y — x J L — y\ x — y y 

with the +- prescription [K(x, y)] + = K(x, y) — S(x — y) J dzK(z, y) and the shorthand notation 

Q{x,y) = sign(l + y)6 1 



1 + yJ \l + y 

Its eigenvalues can be simply calculated for discrete conformal spin n. If one did this for complex 
conformal spin j, one would encounter the same problems as for moments of GPDs. Namely, 
terms proportional to sin(7rj) would appear. The analytic continuation of the discrete eigenvalues 
can, however, be defined without such terms: 



7 7 (0) = C F ( W(j + 2) - #(1) - 3 - - 2 - ) , C F = ^, m = -^lnT(z) . (100) 



This is just the forward anomalous dimensions of twist-two operators, well-known from deep 
inelastic scattering. The evolution of the conformal moments is thus governed by the ordinary 
differential equation 

jxAm^AV) = -^ 7 j%foAV), (101) 

which is easily solved. Equating the renormalization scale \i with the resolution scale Q and 
inserting the solution into the Mellin-Barnes integral (|81|) leads to 

j2 ^ i f c+lo ° , 1 , , f 7 ? (0) f Q2 daa s (a)} , A 2 , . 

^»,,,A»,e>) = 5 jf^ dJ — Pj{Xin)ew ^ _^ij mj( ,, A2 , sS) , (102) 

where the moments mj(j], A 2 , Qjj) belongs to a given GPD at the input scale Q . 
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3.2 Mellin-Barnes representation for amplitudes 



We study next the convolution of a GPD with a given hard-scattering amplitude. To have a 
concrete example at hand we deal here with the so-called Compton form factors that appear in 
the perturbative description of DVCS: 



Hi, a 2 , a' 



p=u,d,s,g 



1 dx 

iT 



2tt 



Cp \Tq) 



0(a 2 s ) 



F p (z,£,A 2 ,/i 2 



.(103) 



For this process the skewness parameter is equal to the Bjorken like scaling variable, i.e., rj = £. 
The same kinematic constraint appears also in the hard exclusive electroproduction of mesons. In 
fact, the result for the Compton factors, we will give below in Eq. ()106|) . can be adopted for this 
process, too. This is trivial to LO and requires some additional work beyond this order. 
Let us consider the partonic Compton form factors to LO accuracy, 



dx 



Ql 



i - x - ie T £ 



Ql 



+ x — ie 



F p (x,£,A 2 ,Q 2 



(104) 



where Q p is the electric charge of the parton. Employing the decomposition © of GPDs into 
quark and anti-quark parts, the partonic Compton form factor (J104)) reads 



A 2 , Q 2 



dx 



Ql 



=F 



Ql 



£ — x — ie £ + x — ie 



[g p (x,e,A 2 ,Q 2 ) + g p (a;,e,A 2 ,Q 2 )] . (105) 



Next we employ the Mellin-Barnes representation ()81)1 for GPDs and perform the momentum 
fraction integration by means of Eq. ()171|) . Then the Compton form factors are expressed in 
terms of the conformal moments 



T 



Ql 

2i 



c+i °° f _! 2^+^(6/2 + j) 

c-ioc r(3/2)r(3 + J ) 



I — 



cos(7rj) =)= 1 
sin(Trj) 



K + m 3 



(e,A 2 ,Q 2 



(106) 



where the sum of mj + rfij is given by the analytic continuation of odd and even conformal 
moments for the vector and axial-vector case, respectively. As we realize by comparing with Eq. 
(|77jl and Eq. the imaginary part of T p is n [q + q] (^, ^, A 2 ) as it must be. The real part of the 
amplitude, given by a principal value integral in the momentum fraction representation, contains 
in the integrand an additional factor tan(7rj/2) and — cot(7rj/2) for the vector and axial-vector 
case, respectively. 

In a conformal subtraction scheme the inclusion of perturbative corrections is straightforward. 
Higher order corrections can be written as convolution of the LO hard-scattering amplitude (J104J) 
with certain kernels, which are conformally covariant. Hence, in analogy to the discussion in 
Sect. 13 .1\ this yields in the Mellin-Barnes representation a multiplication with the corresponding 
eigenvalues, which are known from deep inelastic scattering. For instance, to NLO accuracy the 
partonic form factor 7i p for quarks in the parity even sector of DVCS on a nucleon target reads 
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for \i = Q 



01 
2% 



c+ioo 



^+ l T(5/2 + j) 

r(3/2)r(3 + J ; 



i + tan 



7TJ 



fo + fy] (£,A 2 ,Q 2 ) 
(107) 



with the NLO coefficients 



(i) 



#(i+j)+iUi(j+2)-?+ 



5-2^(j) 

2(i + i)(j + : 



^(i + i) 



l 



+^7 7 (0) {2Si(2j + 3) - + 2) - S l (j + 1)} 
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where the analytic continuation of harmonic sums is defined by derivatives of the T function 



Si(z) = l«z+l)-lKl) 



^(*) = -^-lnr(z). (109) 

The first line in Eq. ()108|) is up to an overall normalization factor the well-known perturbative 
correction to the Wilson coefficients of the structure function F\ in deep inelastic scattering [SB] , 
while the addenda in the second line is induced by the non-forward kinematics. This result 
is verified by a direct rotation from the minimal subtraction scheme to the conformal one and 
coincides with the prediction of the conformal operator product expansion [21^ I59j. 

Another advantage of the Mellin-Barnes representation is that it might be useful for an analytic 
approximation of the Compton factors at smaller values of £, lets say £ < 10~ 2 , which should 
lead to a rather good approximation of the scattering amplitude for the kinematics in Collider 
experiments. Such an approximation has been already studied within the DD formalism, see (HE] 
and references therein, however, it remained restricted to the perturbative LO approximation and 
only the term containing the leading power in £ could be extracted. The main idea here is to shift 
the integration path in Eq. (|107|) to the left, so that one picks up the leading order contribution 
for £ — > 0. Suppose that for the vector case the first singularity on the l.h.s. is a pole at j = ao 
with — l<a < c <l) which might depend on A 2 . We recall that in this case only the analytic 
continuation of odd moments enters and, thus, c can also be chosen to be positive, however, it 
must be smaller than one. The integration path can then be shifted further to the left such that 
— 1 < d < «o an d all other singularities remain to the left of the new integration path, while the 
leading pole contribution (j = ao) is explicitly taken into account. To LO accuracy the partonic 
Compton form factors thus read 
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r(3/2)r(3 + J ) 

This result could even be improved by further shifts of the integration contour. A more systematic 
expansion, requires that also the conformal moments mj(£, A 2 ) are expanded in £ in the vicinity 
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of £ = 0. We remind that £ is related to the Bjorken variable by £ ~ xb/(2 — Xb) and is for 
present fixed target experiments certainly not larger than ~ 0.4. Hence, already the inclusion of 
the 0(£ 2 ) corrections might be sufficient to obtain a good approximation of the Compton form 
factors in this kinematics. 

Unfortunately, the method for a systematic approximation, pointed out above, must maybe be 
refined due to the following complication. Namely, we know that in general the conformal moments 
will have a branch point at £ = and so a Taylor expansion around this point is not possible, after 
the analytic continuation is performed. Even worse, it turns out that the appropriate analytic 
continuation, which guarantees the correctness of the Mellin-Barnes integral, can induce terms 
that lead to a behavior oc rf or even oc rj 2 '-' . Certainly, the latter term requires that we shift the 
integration path to the right, while the former require a consideration of all poles. Corresponding 
to its behavior at j — > oo, we might close the integration path so that the positive or negative axis 
is included. For specific conformal moments it is probably still possible to obtain a systematic 
expansion of the Compton form factors in powers of £ up to a certain order. In general, however, 
this seems to be a serious problem. 

One might naively expect that this issue can be resolved, if one goes back to the definition 
of conformal moments (|3b[): expand first c n (x,£) in £ 2 and take then moments with respect to x, 
which in turn can be simply continued to complex valued j. Such an expansion looks like, see Eq. 



The coefficients c jm vanish for integer value j — n with m > n and so this expansion degenerates 
then into a polynomial 12 jHZ]. in rj of degree n. Hence we can shift from the beginning for each 
individual term the integration path in the Mellin-Barnes integral to the r.h.s., i.e., c — > c + m. 
However, one easily realizes that this procedure will also introduce new poles in the complex j 
plane that are also shifted to the right of the real axis, remaining, however, to the left of the new 
integration path. Finally, it turns out that these poles will remove the power in £ 2 we gained by the 
expansion (jlllj) . With one word, taking the limit £ — > in general conformal moments, leads to 
the correct leading power behavior of the Compton form factors in £, however, the normalization 
might be wrong. 

Let us finally comment on the representation for the scattering amplitude in the case of the 
production of a hadron pair by photon fusion. For the important phenomenological situation that 
one photon is on-shell, it is of course related to the DVCS amplitude by crossing. By means of 

12 Obviously, for non-negative integer n we are dealing with integrals of the type j_-,dx x n ~ m q(x, r], A 2 ), well 
defined for n > m. To avoid divergencies for complex valued j with 5Re j < m — 1, the integral should be defined as 
a contour integral in the complex x plane so that it exist and can be viewed as analytic continuation with respect 
to the variable n — m 
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the crossing relation (|28j) the amplitude follows to LO accuracy from Eq. (jl04|) : 



F P (CW 2 ,Q 2 ) = (1-2()Q 2 





"1 1 


/ dz 


T 1 


/o 


z 1 — z 



%(z,(,W 2 ,Q 2 ). (112) 



Employing the conformal partial wave decomposition (|43|) for GDAs, we immediately find the 
following series for the scattering amplitude 
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As discussed in Sect. 12.21 the conformal moments M n ((,W 2 ) are related to the GPD ones by 
crossing, cf. Eq. (J43j) . Of course, the inclusion of perturbative corrections and evolution effects is 
done in an analogous way as in the case of DVCS. The only difference is that we are now dealing 
with integer n. The sum (|113|) converges, just as the partial wave decomposition for the GDA itself. 
However, the oscillations caused by the polynomials p n (l — 2z, 1) drop out and so its numerical 
approximation is easier to handle. This representation can be directly used in phenomenological 
studies. Knowing an appropriate analytic continuation of M n ((, W 2 ) would also allow to rewrite 
the conformal partial wave expansion ()113|) as Mellin-Barnes integral: 

r rK , W>, <« = (!- 2 C )f f % ^IVl^'VT^ (114 » 

2? J c-ioo r(3/2)r(3 + j) sm(7Tj) 

Here we employed the symmetry relation M n ((,W 2 ) = (— l) n M„(l — (, W 2 ) and that Mj is 
holomorphic in the first and forth quadrant. 



4 Evaluation and parameterization of conformal moments 

This section is devoted to the parameterization of conformal moments and the numerical treat- 
ment of GPDs and Compton form factors. We study first the analytic properties of conformal 
moments for a simple toy GPD ansatz and give then examples for the numerical evaluation. Then 
we introduce a simple parameterization of the conformal moments with respect to the skewness 
dependence, i.e., we consider only a reduced GPD and its crossing analog, and discuss its depen- 
dence on the momentum transfer squared. Finally, we suggest a simple GPD model, which is 
rather flexible in its parameterization. 

4.1 Numerical treatment of the Mellin-Barnes integral 

Let us consider the conformal moments in terms of the Radyushkin ansatz (JT7j) for the reduced 
GPDs. This toy example can be treated for integer values of b and (3 in an analytic manner. To 
represent the results here in an explicit form we choose again 6 = and our toy model with (3 = 1. 



32 



In this case the function u(x, 77) is given in terms of Eq. (jl2J) by 

<K*,l) = t>(-V < * < 1)^ (y±J + % < * < 1)^ (^) ' + ° . (115) 



The analytic continuation of the function /i n (— rj), see Eq. (|38|l . reads 

r(l/2)r(3 + j) n /-j,3+j,2 + a 

= " 2^(3/2 + 3F2 l 2,3 + a 



2, y :il6) 

This function has the behavior needed for the derivation of the Mellin-Barnes integral. However, 
it is not appropriate for an analytic continuation to negative 77, since the argument of the hyper- 
geometric function would then become larger than one, i.e., (—1 + t])/2t] — > (1 +r))/2r), and so its 
value is given at the branch cut of the hypergeometric function. To analyse the situation in more 
detail, we can first express the hypergeometric function in Eq. (|116|) as a sum of two 3F2 functions 
with argument 2t]/ (—1 + 77). This gives 

= (2 + cQ(l-T7) 1+J F (-3- 1, -j, -2- j -a 2 V \ (2 + a)rf^) 
^ V) ' 2(2 + j + a)r l 3 *V -2j-2,-j-l-a 
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21+^(1-^)1+" r(3/2 + j) r(-a)r(3 + j + a)sin(vr[j + a}) ' 

The first term on the r.h.s. contains poles on the real axis at j = 1/2, 3/2, . . . and for a < at 
j = —a, 1 — a, . . . Both of them are cancelled by the second and third term, respectively. 

The continuation of this function to negative values of 77 is obviously not unique, since we have 
branch points at 77 = 0, which leads to phase factors e ±ln U+ 1 + a ) as we ll as e ± l27r u+ 1 ) < Such phase 
factors in fij(t]) would, however, violate the assumptions made for the derivation of the Mellin- 
Barnes representation. To ensure that this does not happen we define the analytic continuation 
by 

J 2ism(nj) 

were the phase factors e compensate the phases which arise from the 77 J terms in the definition 
(J35|) of conformal moments, while taking the "discontinuity" and dividing it by 2ism(wj) ensures 
that no new singularities appear on the real axis and that fij(rj) is bound for j — > 00. In fact the 
properties of Legendre functions imply that the prescription (jll8)l provides the contribution from 
the outer region 

^ - \y ^ (si) i+ * + ^ ,+ " (T Tj m^ < <-> 

where the second term on the r.h.s. arises from a pole at x = rj and cancels the contribution from 
the lower bound of the integral. Hence, we realize that the central region is still ignored in Eq. 
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and its contribution might be taken into account within Eq. (|69|). The term that appears 
there from the pole at x = r\ is already contained in Eq. (|119j) and so only the discontinuity on 
the negative axis contributes. 

Analytic continuation via Eq. ()118|) leads for our toy ansatz to 



(2 + «)(! + „ f-j- 1, — j, -2- j -a 



2(2 + j + a)r ] 3 'V -2j-2,-j-l-a 
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^i+i+a r(l/2)r(l+j) r(3 + a)r(l+j-a)sin(7ra) 
_ 2 i+j-a( 1 + r/ )i+" r(3/2 + j) r(-a)r(3 + j + a)sin(7r[7 + a]) ' 

As before for Hj(—rj) there are no net singularities in the first and forth quadrant of the complex 
j plane, however, individual terms possesses poles on the real axis. The first term on the r.h.s. is 
needed to restore the polynomiality in the sum (i n {v) +A i n°( — v)- However, the third term violates 
polynomiality and must be cancelled by the contribution from the central region, still missing. 
We can restore this missing term from the polynomiality condition in such a way that it is free 
of singularities in the first and forth quadrant of the complex j plane. The sin(7ra)/ sin(7r[j + a}) 
term of the third line, however, generates for j = n = {0, 1, 2, • • •}, a sign alternating series and 
to get rid of it, we separately continue even and odd moments: 

/ , x /x Arv x v 1+j+a r(l/2)r(l+j)r(3 + a)r(l+j-a) , , 

mAnW) = + ,f(-v) + V+Ai + r(3/2 + i) r(- a) r(3 + /+«)) ■ < 121 > 

with a = 1 and a = — 1 in the even and odd sector, respectively. We remark that this result for 
the conformal moments of our toy GPD can be alternatively obtained within the framework 

given in Sect. I2.3T21 Especially, the term we restored from the polynomiality condition arise from 
the discontinuity on the negative real axis in Eq. (|69[) . 

Now let us come to the numerical treatment of the Mellin-Barnes integrals for GPDs (|77|8()j) 
and Compton form factors (jlOfij) . The numerical evaluation can be easily done once the conformal 
moments are known in analytic or numerical form. Corresponding to the behavior of GPDs at 
x 1, they should vanish at large j rather fast [T^]. Hence, the integral converges rather fast, 
too, and so one can in practice perform the integration over a finite interval. If one includes higher 
order corrections or the evolution the numerical treatment remains stable. 

For our toy model we display in Fig. EJa) for a = —1/2 and t] = 0.25 the GPDs arising 
from odd and even conformal moments, while the sum of them (solid line) provides the original 
support. The evolution with respect to the renormalization/factorization scale /i 2 = Q 2 is taken 
into account simply by including the evolution operator in the Mellin-Barnes integral. For the 
vector case the even quark moments evolve separately, i.e., they do not mix with the gluonic ones, 
and so we can employ Eq. ([102)1 . For a s (Q) we take its LO approximation for three quark flavors 
and A = 0.22 GeV. That the numerics is completely unproblematic is demonstrated in Fig. EJb), 
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Figure 6: In panel (a) the toy GPD (|115jl is displayed for a = —1/2 and rj = 0.25 as solid line and 
its symmetric and antisymmetric part as dashed and dotted line, respectively. Panel (b) shows 
the evolution of the symmetric part from the input scale Qq = 0.5 GeV 2 (solid) to the scales 
Q 2 = 1 GeV 2 (dash-dotted), Q 2 = 10 GeV 2 (dashed), Q 2 = 100 GeV 2 (dotted), and (nearly) 
asymptotic limit (thin dotted). 

where the asymptotic case Q 2 — > oo is nearly reached by setting Q 2 = 10 1000 GeV 2 . We remind 
that in this limit the outer region will die out and only the lowest conformal moment contributes 
leading to the asymptotic GPD 

q ^(x, V ) = 6( V -\x\)^-J^^. (122) 

The numerical procedure for the calculation of scattering amplitudes is even easier to handle 
than that for GPDs themselves. As mentioned, a nice feature of the Mellin-Barnes integral is that 
it can be used to derive an expansion in powers of £. Let us suppose that we like to evaluate 
the Compton form factor (J107|) to LO accuracy, where the charge Q p is set to one. Here only the 
analytic continuation of the odd conformal moments is needed. Guided by the small x behavior 
of parton densities, we choose the parameter a to be negative and larger than -3/2. To ensure 
that all singularities are on the l.h.s. of the integration path, we take — Max(l/2, — 1 — a) < c < 
Min(l/2, — a). The first pole which appear in the integrand on the negative axis arises from 
tan(7rj/2) and is at j = —1. We remark that the pole which we would have expected in the 
forward case, namely, at j = — 1 — a is absent for £ > 0. Rather there appears a new one at 
j = — 1 + a, which is associated with the behavior of the GPD at the cross-over point x — rj. 
From our explicit expression for the conformal moments, given in Eqs. (Ill 7)1 . (|12()jl . and (|121)1 it 
is obvious that the integrand contains three different pieces, proportional to and £ Q . 

For the term proportional to £~ 1 ~ J we can arrange a systematic expansion in £ by a shift of the 
integration path to the left. The poles appear here at j = { — 1 —a, —2 — a, ■ ■ •}, j = { — 1, —2, • • •}, 
and at j = { — 1/2, —3/2, ■ ■ •}. In the latter case the contribution will be cancelled by those from 
the piece proportional to which has poles on the positive axis only, at j = {1/2,3/2, • • ■}. 
This is established by a shift of the integration path to the right. What remains are the terms 
proportional to £ a . Here we can close the integration path so that the first and forth quadrant is 
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Figure 7: The Compton form factor H(£,A 2 = 0, <2 2 ), arising from the toy GPD (|115jl . to 
LO accuracy. Left: Exact imaginary (solid) and real (dashed) part of the Compton form factor 
and their approximation by Eq. (jl23j) (dash-dotted and dotted lines, respectively). Right: The 
Compton form factor evaluated numerically at the input scale Qq = 0.5 GeV 2 (solid and dashed) 
and evolved with the flavor non-singlet evolution equation to Q 2 = 4 GeV 2 (dash-dotted and 
dotted). The imaginary part is shown as solid and dash-dotted lines and the absolute value of the 
real part as dashed and dotted lines. 

included and employ then the residue theorem, where the only poles are at j = {—a, 1 — a, ■ ■ ■}. 
All of these poles contribute to the leading power behavior and they must be resummed. Finally, 
collecting the results and neglecting power corrections of order 0(£ 3 ) leads to the approximation 



= 2°ir(2 + a)[i + cot(^)]e 1 {l-i(l+a)^tan(^ 



+ (1 + a)(2 + °V 4- OK 3 )) " 1 1 + , (2 + " + Q V + 0(f) \ . (123) 



a \ (1 — ck)(2 — a) 

This expansion coincides with that for the exact Compton form factor, which for our toy model 
()115j) is exactly calculable. It is remarkable that this expansion does contain odd powers of £. In 
fact, the £ term in the first line should be especially important in the small £ region. Numerically, 
this approximation works quite well for the region that is of phenomenological interest. The 
deviation from the exact expression is about 3% for a = — 1/2 and £ = 0.4 for the imaginary part. 
For the real part the approximation induces a small shift of the zero of the exact expression from 
£ ~ 0.456 to £ ~ 0.431. The accuracy of the approximation ()123j) grows with increasing a and 
rapidly with decreasing £. It is amazing that this approximation remains qualitatively correct as 
long as £ does not approaches one, see left panel in Fig. [7| In the right panel of this figure we 
show the imaginary and real part of the Compton form factor, evaluated numerically with the 
Mellin-Barnes integral, for smaller values of £. Here the difference between the approximate and 
exact result is even invisible. Let us stress specifically, that we do not encounter any numerical 
problems in the small £ region. 
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4.2 Ansatze for conformal moments 



The conformal moments for complex valued conformal spin, appearing in the GPD representation 
(J77|) . can be expanded in terms of 0^(1,77) using a suitable integral transformation. Such an 
integral transformation can be viewed as the analytic continuation of the conformal expansion 
()44)1 for non-negative integer conformal spin. In such an representation the 77 and A 2 dependence 
in the conformal moments separates and only an ansatz for the form factors Fjk(A 2 ) is required. 

Lets have a closer look at these form factors. The GPD Mellin-Barnes integral representation 
(|77l can be interpreted as describing the effective summation of all particle exchanges in the t- 
channel. These are labelled by their conformal spin. Of course, conformal symmetry is broken in 
the non-perturbative sector. However, the conformal spin can still be used for the classification 
of excitations, just like in quantum mechanics for a non-spherical potential, the partial wave 
expansion in terms of spherical harmonics can still be employed to solve the Schrodinger equation. 
For vanishing skewness, i.e., 77 = 0, the quantum number of the conformal spin can be replaced by 
the common spin J = j + 1 and the conformal moments are given by the "diagonal" form factor 



cf. Eq. (|I3J). Having in mind that the conformal partial waves pj(x,r] = 0) are then given by 
(l/x) J , where 1/x plays the role of a (rescaled) energy, the Mellin-Barnes integral for GPDs looks 
similar to the t-channel scattering amplitude at large energies. This suggest that there exists a 
connection with Regge theory or at least with Regge phenomenology. Let us remind that in deep 
inelastic scattering, i.e., for A 2 = 0, this connection shows up in the small x behavior of structure 
functions, which is governed by the intercept of the corresponding leading Regge trajectories 
Note, however, that the small x behavior of parton densities depends on their conventions, i.e., 
on the factorization scheme and scale. A more recent analysis for the unpolarized valence quark 
densities can be found in Ref. [HH]- Below it will be demonstrated that for 77 = there exist 
phenomenological indications that the A 2 dependence in Fjj(A 2 ) is related to Regge trajectories. 

Also for general kinematics, i.e., 77 7^ 0, one can take this duality between t- and s-channel 
serious, see the discussion in Sect. 12.11 A description of DVCS in the high-energy limit is given 
in [70], where the leading Regge trajectory arises from the pomeron exchange in the t-channel. 
Moreover, the t-channel description provides arguments for the so-called D-term [HS] , appearing in 
the nucleon GPDs H (x, 77, A) and E(x, 77, A), and for the so-called pion pole term [7_T], appearing 
in E(x, 77, A). Note that these effects are taken into account by a modification of GPDs that affects 
only the central region and so their s-channel counterpart is absent. On the basis of the conformal 
partial wave expansion, applied to LO accuracy, and the crossing relation between GPDs and 
GDAs a dual description of the former ones in terms of t-channel exchanges has been suggested 
in Ref. |SH1- Here the conformal moments have been decomposed into contributions with definite 
angular momentum. Since the concept of conformal spin has to the best of our knowledge not been 




(124) 
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worked out for applications in hadron spectroscopy, one should focus on this more appropriate 
quantum number. Nevertheless, let us point out that the form factors Fjk(A 2 ) for j ^ k, defined 
in Eq. (JHJ), measure the strength of (non-perturbative) conformal symmetry breaking. In this 
paper, however, we do not proceed with the spectroscopic interpretation of these form factors (or 
some rotated version of them) and the search of an appropriate representation for the conformal 
moments with complex valued conformal spin in terms of them. 

Instead, our aim is now more pragmatic, namely we will introduce and study specific ansatze 
for the conformal moments. Therefore, we now introduce and explore ansatze that are simpler to 
handle. We start in Sect. 14.2.11 with a reduced variable dependence by setting A 2 = 0. Then we 
implement in Sect. !4.2T2l the A 2 dependence. In Sect. 14.2. 31 we especially consider this dependence 
for the case rj = and discuss the resulting GPD H(x, rj = 0, A 2 ), especially, its interpretation as 
three dimensional parton density. 

4.2.1 Ansatze for reduced conformal moments 

For our simple toy GPD example ()115)) . studied in Sect. 14. X| the conformal moments appear to be 
rather complicated functions. Indeed, it is not clear at all whether the complicated structure of 
rrij(r), A 2 = 0) in our toy model, is an artifact of the analytic continuation procedure or indicates 
some physics, related to the skewness dependence. So let us explore several "minimal" ansatze 
for the conformal moments. 

To ensure that the forward limit is correctly reproduced, we factorize the poles in the complex 
j-plane that survive the limit r] — > 0: 

/ -»? i ^ T(a + l+j)T(a + 3 + 2) 1 , 9 , 
m 7 (r7, A = 0\a,/3) = —7 ^J' y n H J -b, (rf) . 125 

Here the normalization of the function is bj(rf = 0) = 1 and, moreover, we considered it as 
convenient to normalize the lowest moment according to m (r], A 2 = 0\a, f3) = 1. In the forward 
limit we arrive after an inverse Mellin transform at 

m ,(, = 0,A* = k/3 ) + (126 ) 

which is well defined for all realistic values of a, except for a = — 1. For this special value one 
can first multiply the expression with an appropriate normalization factor that cancels the factor 
l/r(a + l). The standard parameterization of a parton density, i.e., Nx a (1 + A\fx + Bx) (1 — Xf, 
can be obtained in Mellin space by the linear combination 

qj = N'[ mj (0,0\a,p) + A' m j (0,0|a + 1/2, (3) + B' m j (0,0|a + 1,(3)} , (127) 

where 

r(i + n) r(i + fl r( 3 /2 + a )r(2 + a + f3 ) B , = i + a B (128) 

T(2 + a + p) ' r(l + a)r(5/2 + a + /3) ' 2 + a + (3 K J 
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It is required that the function bj(rf) for non-negative integer values of j — n reduces to a 
polynomial of order n/2 and (n ± l)/2 for even and odd n, respectively. Here the order of the odd 
moments depends on the specific GPD. Moreover, these functions have only singularities in the 
second and third quadrant of the complex j plane. We can allow that these functions grow with 
the real part of j, however, they must be bounded for large imaginary parts of j, |arg(j)| < n/2. 

One can classify the conformal moments b n (rj 2 ) in general and their analytic continuation with 
respect to their dependence on both variables 7? and j, which leads to a classification of reduced 
GPDs. (We do not know whether this mathematical fact is known already from some other 
context.) For the time being, we concentrate on functions bj(rj 2 ) that can be expanded around 
the point rf = 0. Several simple examples can be given in terms of hyper geometric functions 



bj [T] 



{a,b},{r,s},{a,p}) = V 7 — - - r-A (129) 
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( -j/2+(l-2p)(l-<r)/4,a 
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F /-j/2+(l-2p)(l-<r)/4,a 
1 I b 
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where the normalization condition at rj = is satisfied. Here the parameters a and b may depend 
on j. For the analytic continuation of even moments we set as above a = 1 in the case of odd ones 
a = — 1 and the choice p — 1 leads for j = n to polynomials of order (^ 2 )( n+1 )/ 2 ; while p = gives 
polynomials of order (^ 2 )( n ~ 1 )/ 2 . To ensure that no branch cut appears in the interval < rf < 1, 
the parameters r and s must fulfill the inequalities r < 1 and r + s < 1. 

The definition (|129|) is quite general and contains a number of special cases. For instance, if 
we set a = 1, a = (1 — j)/2, 6 = 2, and r — 1, 

hf^Sh -W9 91/1 1 /1 nl^ r(3/2)r(3 + j) f-j/2,(l-j)/2 2 \ 

b j(jl {(1 - j)/2, 2}, {1, s}, {1, 0} j = 2j+ i r ^ 2 + 2Fl ^ 2 ^ J' ^ ' 

we recover for s = — 1 the definition of conformal moments (J35|) with x — 1. Here they are 
expressed in terms of hypergeometric functions that arise from the original ones by a so-called 
quadratic transformation. These moments can be used for even and odd values of j. If one needs 
odd moments that are of order (ry 2 )( n+1 )/ 2 ; one should set a = (—1 — j)/2. We remark that the 
expansion in the vicinity of if = exists only as linear combination of two power series in rj 2 , one 
of them containing the overall factor (t7 2 ) 3 / 2+j . 

Besides the other parameters, s controls the strength of the skewness dependence. If we set it 
to zero, the function bj(rj 2 ) is simply one. 

If a = b, bj(r] 2 ) reduces to the simple function 



bj 77 



/ \(i+P)/2 
{a,a},{r,s},pj = (l -— - V 2 \ . (131) 

Here 1 — r appears as a scaling factor and so in the following r can be set to zero, when simul- 
taneously s is restricted to s < 1. For negative values of s the function (jl29|) is even analytic for 
1 < rj. To get a clue which values of a(j) and b(j) are allowed in the ansatz ()13()j) . we generate 
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cases 


parameters of bj irf {a, 6}, {r, s}, {1, 0} j 


explicit expression 




{a,b} 


{r,s} 




(a) 


{(l-j)/2,2} 


{1,-1} 


r(3/2)r(3+j) p ( -j/2,(1~j)/2 1 2 \ 
2^+^(3/2+^)^1^ 2 ' ) 


(b) 


{a,b} 


{1,0} 


1 


(c) 


{-199/4 + j, 2/3 + 2j} 


{0,-1/4} 


2 ^ 1 ^ 2/3+2j 


-*) 


(d) 


{-199/4 + j, 2/3 + 2j} 


{0,1/4} 


^ 1 ^ 2/3+2J 


*) 


(e) 


{a, a} 


{0,±l/4} 





Table 1: The parameters and resulting functions for our ansatze concerning the analytic contin- 
uation of even conformal moments f)129|) . classified by 



"associated" conformal moments of bj iff {a, a}, {0, s},p) by the convolution integral 



dzfj(z)(l - zsr] 2 ) {j+p)/ \ with dzfj{z) = l. 



(132) 



Here it is required that fj(z) as function of j is bound for j ^ oo with |arg(j)| < tt/2 for < z < 1. 
To arrive at the parameterization in terms of hypergeometric functions we choose 



/;(*) 



T(a(j))T(b(j) - a(j)) 



(133) 



The requirement for the bound of fj{z) is certainly satisfied for 5te a(j) > and 3?e(6(j)— a(J)) > 0. 
Especially, if a(j) tends to infinity for j — * oo, b(j) has to grow as a(j) or even faster. The case 
a = (1 — j)/2 and b = 2, mentioned above, can be obtained within an analogous treatment. Here, 
however, one has to choose an integration contour in Eq. (|132|) in the complex plane that encircles 
the point z = and so the convergence condition for the integral can be relaxed. We will skip 
this issue here and refer, for instance, to Ref. |d7] . 

We now explore the resulting GPDs and GDAs. As example we take the reduced valence quark 
GPDs and GDAs in the vector case with the rather realistic values a = — 1/2 and /3 — 3. The 
unpolarized parton density is normalized to one and reads, cf. Eqs. (J125J) and 



qix) 



32 v 



x) 



(134) 



The function bj{rf) is the analytic continuation of conformal moments with even n and the pa- 
rameters are specified in Tabled A few comments are in order. In case (a) we took the conformal 
moments itself, here in a more appropriate representation that is for rj > equivalent to Cj(l,rj). 
They possess the properties we required in the derivation of the Mellin-Barnes integral and of 
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course, they reduce to polynomials for both even and odd non-integer values of n. We will use them 
in the Mellin-Barnes integral (J79|) within even conformal partial waves (J75j) . Crossing symmetry , 
see Eq. ()4llj1 . requires the existence of an analytic continuation to 77 > 1 and this is achieved here 
by the change of arguments Cj(l, 77) — > Cj(l — 2(, 1), see Eq. (|44|). For the numerical evaluation of 
the GDAs the conformal partial wave series f)43j) or alternatively the Mellin-Barnes integral ()84|) 
can be used. The projection on the even moments can be achieved by symmetrization: 

$(z,C)-i[$(z,C) + $(l-z,C)] • (135) 

The case (b) seems to be trivial. But after crossing it results in (1 — 2()i , which leads for odd 
moments with ( > 1/2 to an alternating series with an absolute value that converges to one in the 
end-points. For the convergence of the conformal partial wave series f)44j) in terms of polynomials 
an exponential suppression factor 2~ n is required. Hence, this series only converges in the interval 
l/4<C<3/4. The analytic continuation with respect to ( can be achieved by the Mellin-Barnes 
integral ()85|) . However, outside of the convergence region, i.e., £ < 1/4 or 3/4 < (, we will find a 
GDA that does not vanish at the end-points z = {0, 1}. 

The cases (c) and (d) differ by the sign of the argument in the hypergeometric function and 
follows from (e) by an integral transformation (|132j) . They reduce to polynomials for even n only. 
The factor 1/4 in the argument improves the convergence property after crossing. The second 
parameter in the upper line of the hypergeometric function, i.e., — 199/4 + j, can for certain non- 
integer values of j = n < 49 be a negative integer. As explained above, this does not generate 
problems. The choice of this "big" constant —199/4 induces a numerical enhancement of 0(rf) 
terms. The parameter 3/2 + 2 j in the lower line compensates for rather large values of j the growing 
of the second argument in the first line. We did not include any poles on the positive j axis. As 
a consequence, we have new poles on the negative one. They appear at j = —1/3, —4/3, • • ■ and 
die out in the limit 77 — > 0. With our choice a = —1/2, which determines the small x behavior of 
the parton densities and is associated with a pole at j = —1/2, a new "leading" pole at j = —1/3 
arises for non-zero skewness. Its contribution, however, should be suppressed as 0(rf). So for 
instance, in the Compton form factors it should only give rise to a r] 2 ri~ 2 ^ = 77 4//3 term. Moreover, 
we changed the normalization of the other poles for < 77. For not too large values of 77 this 
should produce only a numerically small effect. Since we included a numerical enhancement, this 
might induce some sizeable changes of the GPD and GDA shapes for not too small values of 77. 

In Fig. [Sj we display the 77 and ( dependence for the second moment, i.e, of 62 iv 2 )- The 77 
dependence of the conformal moment in case (a) [(solid line)], i.e., 

czM) = 1-^ => c 2 (l-2C,l) = (l-2C) 2 -^, (136) 

is rater weak and so it is "similar" as in case (b) [constant or (1 — 2Q 2 (dash-dotted line)]. In 
contrast we find for cases (c) [dashed line] and (d) [dotted line] a large deviation in opposite 
directions. For GPD moments it will die out for 77 — > 0, while for the GDA moments the difference 
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Figure 8: The r\ (left) and ( (right) dependence of b n {rf) and (1 — 2Q 2 b n {l/{l — 2() 2 ), respectively, 
for n = 2 are shown for the parameterizations given in Table ^ (a) solid, (b) dash-dotted, (c) 
dashed, and (d) dotted lines. 




Figure 9: The momentum fraction shape of reduced GPDs with 77 = 0.5 (left) and GDAs with 
( = 0.4 (right), same labelling as in Fig. |H| 

is nearly ( independent and large. We will skip here the detailed discussion of the dependence 
in higher moments, which are suppressed by the factor B(l/2 + n, 4 + n) ~ 1/n 4 for n — * 00. 
Note, however, that only the conformal moments (a) for rf > continuously tend to zero with 
increasing n. Also the rescaled conformal moments (a) for all values of ( possess this behaviour. 
We remark that the conformal moments (e) will mostly not give quantitatively different results as 
(b), so we will in the following not present them. 

In Fig. El we show the resulting GPDs for 77 = 0.5 and GDAs for ( = 0.4. As has already 
been discussed the skewness dependence of the conformal moments in Fig. |H1 cases (a) and (b) 
can be hardly distinguished, for GPDs they are quite the same. In case (c) the area in the central 
region and the magnitude of the maximum at 77 = are enhanced, while the area of the outer 
region shrinks (the lowest moment of all functions is normalized to one). Consequently, also the 
magnitude at the cross-over points x = ±77 decreases. The reverse situation is observed for case 
(d). This is caused by the analytic properties of the conformal moments with respect to the 
variable j, which we explained above. Comparing both panels in Fig. El one immediately sees 
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Figure 10: The reduced GPD at the cross-over point x — r], multiplied with the factor ir, (left) 
versus r] and the convolution of the GDA with 1/z (right) versus (, same as in Fig. (JSJ). 

that the "width" of GPDs shows up in the end-point behavior of GDAs. It is remarkable that 
the GDAs (a) and (b) possess shapes that are governed mainly by the first partial wave of the 
conformal expansion. Indeed, for case (a) it turns out that in the region 0.05 < ( < 0.95 the higher 
partial waves give at most a 10% percent correction. In the much narrower region 0.3 < ( < 0.7 
this is within the same accuracy also true for the GDA (b). This region starts slightly about or 
below the value, where the series turns over to be divergent, i.e., for ( < 0.25 and 0.75 < (. A 
slightly smaller convergency radius holds for the GDA (c) and (d). 

In Fig. El we show the quantities that enter the scattering amplitude for a hard exclusive pro- 
cess, in which only the virtuality of the incoming photon can be varied. Under such circumstances 
the Compton form factors, i.e., Eq. ()104|) or their crossing analogues Eq. ()112|) . are measurable. 
The non-perturbative distributions can not be directly determined by deconvolution, but at least 
within our toy ansatze the correspondence between momentum fraction dependence and the ampli- 
tudes is clear-cut. In the left panel of Fig. E3we depict the GPD at the point x = rj multiplied by 
7T. This cross-over point trajectory is accessible in single spin asymmetry measurements. Within 
our ansatze both the normalization and the slope encodes information about the GPD shape. 
Compared with the left panel in Fig. one realizes that the narrowest (widest) GPD has the 
smallest (largest) value at the cross over point x = rj. The differences of the trajectories diminish 
with decreasing rj. Cases (a) and (b) are only distinguishable when t] tends to one. Certainly, the 
slope of the cross-over point trajectories is dictated by the strength of the skewness dependence 
of the conformal moments. 

In the right panel of Fig. El we show the ( dependence of the scattering amplitude, given as 
convolution of a GDA with a LO hard-scattering amplitude, see Eq. ()112|) . It has been evaluated 
within the partial wave expansion ()1 13j) . taking the first thirty terms into account. The result is 
only displayed for the region in which also the series for the conformal moments (b)-(d) converges, 
i.e., for 0.25 < £ < 0.75, The end-point behavior of GDAs determines the normalization of the 
scattering amplitude while the ( dependence is in all cases (rather) flat. As already mentioned, 
these convergency problems arise from an exponential growth caused by the factor 2 n that is 
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contained in the normalized conformal partial waves (|39|). One would expect that the conformal 
partial wave expansion of GDAs converges for all physical values of ( and z. To ensure this, the 
conformal moments must be exponentially suppressed by a factor 2~ n for large n. We remind that 
the normalization has been adopted from parton densities and so this factor will drop out in the 
limit 7] — > 0. 

We come now to the t] dependence of the cross-over point trajectory for smaller values of r]. 
The trajectory is represented by the Mellin-Barnes integral 

1 r + '°° J ,..-,-i 2,+ ' r ( 5 / 2 + j) r(i/2 + j)r(9/2), 2 



,( "' "» = a t ra r ( 3/2)r( 3 + J ) wf5+ ^ > ' (137) 

which can be used to derive, as explained in Sect. 13 .2\ an analytic expansion in terms of powers 
in rj. The leading terms 

are "universal" in all cases and arise from the poles at j = —1/2 and j = —3/2. The accuracy 
of this approximation for the cross-over point trajectory is for all cases about 1% for rj = 0.05 
and, of course, starts to be much better with decreasing rj. Corrections to the approximation 
(J138J) depend on the specific model and will not be discussed here. For large values of rj they are 
smaller for ansatz (a) than for (c) or (d). Remarkably, in case (b) the integral (|137|) can be exactly 
calculated: 

, , (2-T7) 7 / 2 A/2, 5/2 2\ , . 

Let us summarize the lessons from this investigation. Unluckily, after crossing it turns out 
that three of our models, (b)-(d), lead to convergence problems for the conformal partial waves 
series. Since these series should converges, the conformal moments must for r] > 1 exponentially 
decrease as 2~ n with increasing n. Viable ansatze for the conformal moments are holomorphic 
functions of the conformal spin j that respect the following requirements: 

• they should be bound for j — >• oo with |arg(j)| < 7r/2 for \rj\ < 1 

• An expansion of these functions in rj 2 should exist 

• they should show an exponential suppression for rj 2 > 1 by a factor 2~ n (for integer n) 

In addition they should be flexible enough to be able to generate a large variety of GPD shapes. 
In our examples only case (a), i.e., the analytic continuation of Gegenbauer polynomials, satisfy 
the three requirements. Applying suitable integral transformations to them allows to generate 
"associated" polynomials. We did not study in detail how flexible the GPD shapes can be param- 
eterized within this method but got the impression that only small changes for large £ are possible. 
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We are award that these examples do not cover all possible types of conformal moments. One 
important lesson is that the analytic properties of conformal moments determine the qualitative 
features of GPDs. Certainly, our examples here give only a first insight into this connection. For 
completeness, we mention the existence of a symmetry relation that arises from the definition 
(|35|). namely, the invariance under the replacement j — > — j — 3 



2 1+j r 3/2 + j . . 2 1+ -T3/2 + j _, . . n An . 

, , ', -rn 3 cAx,n) = - - \ ' J -4- ri J cAx,ri) . (140) 

r(3/2)r(3 + j) 3K ' /; r(3/2)r(3 + j) j->-j-3 v 1 

Since, however, perturbative QCD, e.g., anomalous dimensions, already violate this symmetry, 
there is no need to implement it in a non-perturbative ansatz for conformal moments. 

4.2.2 Implementation of momentum squared dependence 

Now we address the implementation of the A 2 dependence. We might include it by introducing A 2 
depended parameters in the ansatz ()125j) in such a way that no singularities in the first and forth 
quadrant of the complex j— plane appear. For the sake of simplicity we write here a factorized 
ansatz 

, a2i n\ T(a + 1 + j)T(a + f3 + 2) , 

™An,A>,P) = r(o + i)r(o + <? + 2 + j) t,tf)F,(A) ' (141) 

with the normalization condition Fj(A 2 = 0) = 1. In the case that the form factors Fj(A 2 ) are 
j independent the x and A 2 dependence is obviously factorized, too. Lattice calculations for the 
first moments of u and d quark GPDs suggest that such a factorization is actually not correct. 
According to these results, the cut-off mass squared in a dipole fit increases with j = n. Unluckily, 
the systematic theoretical uncertainties of these results, especially, those associated with the chiral 
extrapolation are still kind of large, but the preferred fits suggest a linear growth with (j + 1) 13 
|72j . A second constraint on the A 2 dependence arises from the lowest moment of GPDs, which 
is related to partonic form factors, e.g., for the GPD H: 

Ff val (A 2 )= f dxH Pml (x, V ,A 2 ). (142) 

Where Ff' val (A 2 ) can be expressed in terms of the proton (p) and neutron (n) electromagnetic 
Dirac form factors according to 

2F 1 u ' val (A 2 ) = 2Ff(A 2 ) + 2F x n (A 2 ) , F^A 2 ) = Ff(A 2 ) + 2F 1 "(A 2 ) . (143) 



(n + \ )m 2 d q dependence for the squared dipole masses m 2 d , taken from Table I of Ref. [72] gives a 
very good fit for the generalized form factors A u n ~ Q d (m n = 897 MeV ) and A u n ^ d (m^ = 744 MeV ) with n = 0, 1, 2. 
For the large pion mass = 897 MeV one might have the impression that the growth is stronger (weaker) for 
the (axial-)vector case, which might be caused by a different intercept in the power behavior. Certainly, a definite 
conclusion can not be drawn from present lattice measurements. 
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We remark that the lowest moment of a GPD is simply obtained by setting j = in the ansatz 

dun). 

Now we are in the position to build a minimal GPD ansatz for the valence quark GPDs H 
which satisfy the theoretical constraints in the forward case and provide for the lowest moment the 
correct A 2 dependence. Suppose we use the parameterization of the forward parton distribution 
at a given input scale in the form (j!27jl . the conformal GPD moments read then, e.g., for the u 
quark 

m^( V ,A 2 ) = J [m j ( V ,A 2 \a,P) + A'm j (r 1 ,A 2 \a + l/2,P) + B , Tn j (r ] ,A 2 \a+l,P)] , 



T(a + l + j)T(a + (3 + 2) u 

2\ T-iW.val/ a 2 



A |a - « = ^rm^m^ ^ (A/(J + [)) • (144) 

normalized to 2 for j = and A 2 = 0. Guided by the lattice results, we rescaled here the dipole 
masses in the form factor i^"' val by j + 1, which results in a dependence on the ratio A 2 /{j + 1). 
The remaining free parameters A',B',a,(3 are taken from the parton density parameterization 
and for bj{rf) one can use one of the functions suggested above, where preference is given to case 
(a). Analogously, one can deal with other species of GPDs that reduce in the forward limit to 
parton densities. However, for these the A 2 dependence is typically far less known. In particular, 
unpolarized sea quark GPDs are not constraint by elastic form factor measurements. 

Let us have a closer look at the parameterization (|144j) . For j = we obviously find from this 
parameterization 

f dxH Uv Jx, r), A 2 ) = 26 (r / 2 )F 1 "' val (A 2 ) , (145) 

J-ri 

where the ^-dependence drops out, bo(r] 2 ) = 1 and the correct normalization is ensured. Setting 
A 2 = 0, it is also clear from Eq. (jl41|) that we arrive by an inverse Mellin transform at the parton 
density. It has been mentioned above that the GPD for 77 — > and small x is dominated by the 
leading Regge trajectory. In the following we argue that after a small modification this trajectory is 
already present in our parameterization. In our case we are dealing with the p° and u trajectories, 
which according to the analysis of Ref. are parameterized by a linear A 2 -dependence 



a^(A 2 ) = 0.42 + A 2 0.95GeV" 2 , a p (A 2 ) = 0.48 + A 2 0.88 GeV" 2 . (146) 

If we take for the elastic nucleon form factor the standard dipole parameterization and consider 
only the first pole in the complex j-plane of the parton density Mellin-moments, i.e., the pole 
that appears on the negative axis at the largest value of j, our parameterization for the conformal 
moments reads 

A ') « \—. 2T- ?' + 1 + aL T ( 14T ) 

+ J m d(j + 1 + a) — A 2 (l + a) 



3 a 



A B 
+ 



4M 2 N - (1 + c)A 2 /(j + 1 + c) m 2 - (1 + d)A 2 /(j + 1 + d) 
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Here is the nucleoli mass, rrid is the dipole mass appearing the parameterization of the Sachs 
form factors, e.g., in the electric one of the proton G P E (A 2 ) = 1/(1 — A 2 /m^) 2 . A and B are two 
constants, the values of which can be read off, e.g., from Ref. Moreover, we modified here the 
rescaling of the A 2 dependence in such a way that the pole at j + 1 + a = cancels against the 
scaling factor (j + l + a). The constants c and d, appearing in the scaling factors of the remaining 
A 2 dependence, are chosen to be positive so that they do not interfere with the leading Regge 
trajectory. This linear trajectory is given by 

1 -I- n 

a(A 2 ) = a(0) + a'(0)A 2 , a(0) = -a , a'(0) = (148) 



and so we find with Eqs. ()144j) and (j!47j) 

vai, ^ = 2 T(a + 2 + j)T( a + P + 2) /^(A 2 )^ 2 ) 
j KV ' ' 1 + A' + B'T(a + l)T(a + f3 + 2 + j) (j + 1) - a(A 2 ) { 



(149) 



where the remaining A 2 dependence is accumulated in the function 
2 A B 

^' (A ) = AM 2 N -(l + c)A 2 /(j + 1 + c) + m 2 - (1 + dWTO + 1 + d) ' (150) 

As we saw, after rescaling of the dipole masses the leading pole at j = — 1 — a in Eq. ()144l) is 
in the parameterization ()149j) shifted to j = — 1 + a (A 2 ) by an amount proportional to A 2 . The 
numerical value of the dipole mass is rrid = 840 MeV and that of a depends on the factorization 
scale. When we choose the scale to be the intercept of the Regge trajectories ()146|) we find the 
following slopes 

a4(0) = 0.82 GeV" 2 , a' p (0) = 0.74 GeV" 2 . (151) 

These values are only about 15% smaller than the ones given in Eq. (|146|) . In view of the fact 
that we used just the standard dipole parameterization of the elastic form factors and neglected 
all non-leading Regge trajectories, this agreement is quite astonishing. We interpret it as further 
evidence that Regge theory or at least Regge phenomenology is indeed applicable to GPDs. A 
deeper understanding of this issue could provide the key to a dual description of GPDs in terms 
of hadronic degrees of freedom and certainly warrants dedicated efforts. 

4.2.3 Numerical consequences for the probabilistic interpretation of GPD H 

As mentioned in the introduction, GPDs possesses for rj = a probabilistic interpretation in the 
infinite momentum frame. In particular the Fourier transform 14 

H g (x, b) = y ^ e~ fb K H q (x, r, = 0, A 2 = -A 2 ) (152) 



14 To avoid a confusion with the definition of the GPD H q (x, rj, A 2 ) in the axial-vector case, we omit the tilde 
symbol for the Fourier transform of H q (x, rj, — A 2 ). The quantities in the two-dimensional impact space b arc 
indicated by their argument b. 
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Figure 11: The w-valence quark GPD H u (x, r] — 0, A 2 , Q 2 ) of the proton is displayed in panel (a) 
versus x for fixed A 2 = —0.25 GeV 2 with a A 2 -dependence arising from a Regge trajectory via 
Eqs. (I148J) . f)149|) . and (|150|) (solid and dash-dotted) and with a rescaled l/(j + 1) A 2 -dependence 
via Eq. ()144j) (dashed and dotted). The GPD is evolved from the input scale Q 2 = 0.5 GeV 2 (solid 
and dashed) to Q 2 = 10 GeV 2 (dash-dotted and dotted). In panel (b) we show for the A 2 / (j + 1)- 
depended GPD H u the expectation value for the square of the impact parameter (j 154)1 . The 
different lines correspond to Q 2 = 0.5 GeV 2 (solid), Q 2 = 1 GeV 2 (dash-dotted), Q 2 = 10 GeV 2 
(dashed), Q 2 = 100 GeV 2 (dotted) and Q 2 = 10 100 GeV 2 (thin dotted). We used the following 
parameter set: A' — B' — 0, a — — 1/2, /3 = 3 and c = d = 1. 

is the probability to find a quark species q inside the nucleon with momentum fraction x at 
impact parameter b. The latter is defined relative to the center of momentum of the hadron, i.e., 
b is the distance of the active parton in the transversal direction from this center. It is worth 
mentioning that the definition of such a center is based on the existence of a Galilean subgroup of 
transverse boosts in the infinite moment frame [73] . Within light-cone quantization such transverse 
boosts have a field theoretical definition in terms of two conserved charges, expressed by the plus 
components of the energy momentum tensor. Their eigenvalues are good quantum number that 
labels the states of the hadron. Conveniently, they are chosen to be zero for the center. 

This probabilistic interpretation of the GPD (j 152)1 has inspired several authors to build GPD 
ansatze to get a first glimpse of the three-dimensional tomographic picture of the nucleon. This has 
often been done using an exponential ansatz for the A 2 -dependence, which might serve it purpose 
in the space like region, however, violates the analytic properties of scattering amplitudes etc. such 
that crossing relations become meaningless. In the following we study our conformal moments 
for t] = under this aspect. The only uncertainty which is left in our GPD representation is the 
j-dependence of the form factor. In Fig. (jll)) we present the momentum fraction dependence for 
fixed A 2 = —0.25 GeV 2 and two different resolution scales. Obviously, the implementation of the 
Regge trajectory (solid and dash-dotted lines) results in a flatter x-dependence compared to a 
simple rescaling of the dipole mass squared with j + 1 (dashed and dotted lines). 

As a side remark, we comment on the factorized A 2 ansatz for GPDs. Although, it is wrong 
in principle, this does not necessarily imply that all estimates for observables fail completely, at 
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least not for smaller values of A 2 . For the ansatze of conformal moments, we discussed in the both 
previous sections, the normalization of the resulting GPDs between the different versions of A 2 - 
dependence varies for 0.08 < x < 0.2 and |A 2 | < 0.3 GeV 2 , i.e., in the fixed target kinematics, not 
larger than 20%. Certainly, for lower or larger values of x the differences especially in the overall 
size can increase drastically, compare, e.g., the solid and dashed line in Fig. ^2( a )- On the other 
hand the suppression introduced by the Regge motivated ansatz (see solid and dash-dotted lines) 
is welcome to suppress sea quark and gluonic contributions, which notoriously are overestimated 
in the factorized A 2 ansatz for hard exclusive electroproduction processes. 

The average distance from an active parton to the center of the nucleon is defined as jHJ 

^ ^ - fSxt - ^ ^ " - °- A " Q2) l- ' (153) 

Within the parameterization (f 144J) . where for simplicity we again rescale the A 2 dependence by 
+ 1), this average distance can be exactly calculated and expressed in terms of forward parton 
distributions 

(g%(s, Q 2 ) = ) L 4-^lnFf(A 2 ) . (154) 
q(x, Qr) oA z a 2 =o 

The l/(j + l) factor that arises in Mellin space from the differentiation with respect to A 2 gives in x 
space rise to the integral J^dy/y ■ ■ ■. Obviously, by a more refined rescaling of the A 2 dependence, 
see Eq. (j!47|) . we can express the resulting average by a more complex integral. For small x this 
quantity tends to a constant, depending only on the resolution scale, while at large x it vanishes 
as (1 — x). The latter behavior is a consequence of the linear j'-dependence of the dipole masses. 
Such a behavior has been rejected in Ref. [201 - Namely, the quantity 



is interpreted, based on a partonic picture, as the distance of the active parton from the center 
of momentum of the spectators and should therefore be finite for x 1. Consequently, (b 2 ) q (x) 
should vanishes at least with (1 — x) 2 for x — > 1. Although this argumentation is further supported 
by perturbative QCD arguments jTH], it is in our opinion not excluded that the simple parton 
picture can be misleading for quantities, which have no well-defined field-theoretical analog. Note 
also that a (1 — a;) 2 behavior of (b 2 ) q (x) for x — > 1 requires a (1 + j) 2 growth of the squared dipole 
masses. This means in terms of Regge phenomenology that the trajectories should possesses a 
(small) non-linear term. 

Let us come back to our ansatz. To include the scale dependence we use for the GPD H q (x, rj = 
0, A 2 , Q 2 ) in Eq. (J153|) the Mellin-Barnes integral and insert the evolution operator, compare with 
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Eq. dun, 



(0) 

7 j ft 
Jc 



f c+ioo di x ~ j ~ lT 0-+i+ a ) e I ' " u o ° j r, 
(P) q ( x , Q 2 ) = Jc ~ lo ° ni+j+a+^a+i) _ _4_^_ lnjF ^-( A 2) o . (156) 
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a;-j-ir(l+j+a) J 2 ^Sg 
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7 J f Q 2 dcr as(g) 
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For experimental accessible values of Q 2 and for small x, both the numerator and denominator 
are dominated by the leading pole at j = — 1 — a. Shifting the integration path to the left we find 
for x — > the constant value 

= o, Q 2 ) = -^ lnF i val ( A2 )L - ( 157 ) 

We remark that a Regge motivated ansatz would induce a double pole in the numerator and, 
consequently, a logarithmic modification with respect to both the x- and Q 2 -dependencies. Within 
the dipole ansatz for the partonic form factors, taken from Ref . [HE] , we have 

(b 2 } u (x = 0, Q 2 ) = --0.4fm 2 (b 2 ) d (x = 0, Q 2 ) = --0.53fm 2 (158) 

a a 

for the u and d valence quarks. From Fig. ITlT b). where we used the generic value a = —1/2, we 
can read off the qualitative x-dependence of (b 2 ) qml (x, Q 2 ) for given resolution scale Q 2 . It can 
be roughly approximated by a logarithmic growth with decreasing x which changes slope at some 
"cusp" point x cusp (<2 2 ) ~ 10~ 3 . With increasing Q 2 this cusp is washed out. The increase with 
ln(l/x) visible in Fig. II II (b). will not continue in the limit x —> but will saturate at a finite value 
(|158|) . Numerically, we checked this qualitative feature up to Q 2 = 10 10000 GeV 2 and 1/x = 10 40 . 

So within our simple ansatz, which can be easily refined, we reproduce following well-known 
grand picture. At large x the valence quarks are in the center of the nucleon. With decreasing x 
they become ever more delocalized in the transverse direction and will move away in transversal 
direction with decreasing momentum fraction x. In the valence quark region with x ~ 0.3 the 
average transverse distance is of the order ~ 0.4 fm. With increasing resolution scale this value is 
slowly decreasing. In the small x region their transverse distance grows, reaches (for the ansatz 
we used) ~ 0.9 fm in the limit x — > for u valence quarks and a slightly higher value for d quarks. 
Hence, the charge of the proton, which is probed in a scattering experiment, has a non-trivial 
transverse distribution. 



5 Summary and conclusions 

In this paper we have derived a new representation for leading twist-two GPDs and GDAs in 
terms of Mellin-Barnes integrals over the complex conformal spin. This representation is rather 
analogous to the partial wave expansion of scattering amplitudes with respect to complex angular 
momentum, given in terms of Legendre polynomials or functions, respectively. Indeed, also our 
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partial waves can be expressed in terms of associated Legendre functions of the first and second 
kind. Mathematically, there should exist a one-to-one relation to other representations that are 
based on conformal symmetry. The transformations between different representations, however, 
still have to be worked out in detail. For instance, taking the Mellin transform of the Mellin-Barnes 
integral must lead to the integral kernel that maps the "effective forward parton distribution" to 
GPDs |36]l37] or employing the Fourier transform one must arrive at the light-cone position space 
representation [HB1 E3J EH] • We confirm and generalize the results for the gluonic sector given in 
Ref . [3H] . The advantage of our new representation is that the central and outer region are obtained 
from the same conformal moments of a GPD and polynomiality is manifestly implemented. This 
has not been done so far in the approaches, we mentioned. Sometimes, the central region was even 
treated incorrectly, for comments see Ref. jHJJ. We derived a "spectral" representation of conformal 
kernels with complex valued conformal spin. However, we are not completely satisfied with this 
representation, since some support restriction have to be fixed explicitly by step-functions. So far 
an orthogonal eigenfunction basis for these kernels can only be given as series of mathematical 
distributions, which are labelled by non-negative integer conformal spin. 

To leading order the kernels and hard-scattering amplitudes respect conformal symmetry and 
so we could provide the solution of the evolution equation and the Compton form factors for DVCS 
as Mellin-Barnes integral. This representation allows a simple and stable numerical evaluation 
of these quantities. This has some practical advantages, especially, having efficient numerical 
routines at hand for the evaluation of the Mellin-Barnes integral, one might be able to evaluate 
this quantities in "real time" rather than using a database of Compton form factors or GPDs. We 
even demonstrated for the Compton form factors that by means of the Mellin-Barnes integral a 
systematic analytic approximation in powers of £ is feasible. 

Beyond LO order the conformal symmetry is broken in a subtle way by the minimal subtraction 
scheme, applied to the divergencies of composite operators, and the trace anomaly of the energy 
momentum tensor. The first effect can be cured by a finite renormalization, while the latter one can 
be absorbed into either the hard-scattering amplitudes or the evolution of GPDs. Such symmetry 
breaking effects have been studied for 77 = 1 in connection with the pion distribution amplitude 
and the pion elastic and 77* — > it transition form factor. Depending on the model for the pion 
distribution amplitude one finds a 10% variation of the NLO corrections [73]. For skewness 77 < 1 
one would expect an even smaller effect because the conformal symmetry breaking is suppressed by 
a factor rf. For realistic experiments typically 77 < 0.4. Further studies are desirable to clarify this 
issue and, hopefully might provide a systematic expansion of the conformal symmetry breaking 
effects in powers of rf . 

Certainly, the evaluation of the analytic continued conformal moments from a given GPD 
ansatz is a rather non-trivial task, the difficulty of which depends on the analytic properties of 
the corresponding GPD. Since, however, GPDs are almost unknown non-perturbative functions, 
one should rather model the conformal moments directly instead of the GPDs or DDs. This is 
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a non-trivial task with respect to skewness dependence, since certain analytic properties of the 
conformal moments must be respected. Some effort is required to tune the internal skewness 
dependence, arising from the conformal partial waves, accordingly. In the end, however, it turned 
out that our ansatze generated GDAs with non-vanishing end-point contributions. This illustrates 
again that known constraints for GPDs are very difficult to fulfill. An important observation one 
can draw from our examples is that "strong skewness" effects already show up in the first few 
conformal moments for non-negative integer conformal spin. So it is worthwhile to calculate them 
on the lattice to get a clue for the strength of the skewness dependence of GPDs. 

A partial wave expansion of the conformal moments itself avoids the problems described above 
yields a most flexible parameterization of GPDs in terms of form factors. These form factors are 
related to particle exchanges in the t-channel, but this connection still has to be worked out in 
detail. For rj = there is no difference between conformal spin and ordinary spin and the conformal 
partial wave expansion turns over into one with respect to spin. In this kinematical domain we 
found some evidence that Regge phenomenology can be used as a reliable guide for modelling 
the conformal moments. Especially, for unpolarized valence quark GPDs the parameterization 
of parton densities and elastic electromagnetic form factors can be unified and interpreted as a 
leading Regge trajectory. Since GPDs depend on the factorization conventions, such a connection 
can only hold approximatively. Moreover, Regge theory is only applicable for physical amplitudes 
and going beyond the leading trajectory has its own difficulties. Nevertheless, the Regge analogy 
provides some guidance for the modelling of conformal moments. This is especially important for 
those GPDs that are difficult to extract from experiments. 

The advantage of the Mellin-Barnes representation has been demonstrated for several analytic 
and numerical examples, especially, for unpolarized valence quarks in the rj = case. The only 
unknown in this limit is the spin, i.e., j, dependence of the form factors, which for j = 0, are 
measured in elastic electron proton scattering. The momentum fraction dependence of the GPD 
follows then from the t- and Q 2 -dependencies, where the boundary condition at the input scale 
Qq and at A 2 = can be simply taken from the parameterization of Mellin moments for parton 
densities. Remarkably, no additional fitting procedure is needed to satisfy the GPD constraints. 
However, different parameterization of the form factors with respect to j, can lead to a different 
holographic picture of the nucleon. Certainly, the important task here is to pin down the remaining 
degrees of freedom for the j dependence. No question, improved lattice calculation with a realistic 
pion mass can provide at least a partial answer. 

We would like to add a speculation concerning the experimental access to this dependence 
in hard exclusive reactions. Suppose it turns out from lattice measurements that the skewness 
dependence of the conformal moments is weak, the ^-dependence for the scattering amplitude is 
(approximately) known and the only degree of freedom left is the unknown j-dependence of these 
form factors, which determines the shape of the trajectory of the cross-over point of a GPD as 
function of £, in dependence of A 2 . Such a trajectory can be explored in single beam spin experi- 
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ments. Taking Mellin moments of this trajectory, which requires of course some interpolation, one 
directly gets the form factors, up to some normalization factor, in dependence of the (conformal) 
spin. 

Let us finally stress that the crossing relation between GPDs and GDAs is very simple for 
conformal moments. Apart from a trivial rescaling procedure of the skewness dependence it 
involves only the A 2 or W 2 dependence. One has to perform an analytic continuation of the form 
factors from the space- to the time-like region, which requires only a suitable parameterization 
in terms of rational functions (linear combination of monopole or dipole forms), which scale for 
A 2 — > — oo such as predicted by dimensional counting rules. As there exist also non-perturbative 
scales, like the hadron masses themselves or Aqcd, anomalous, i.e., logarithmical deviations, from 
the canonical scaling should be present to some extend. Although, we were not able to give for 
all examples of conformal moments an unified representation of GPDs and GDAs in terms of a 
Mellin-Barnes integral this is not a restriction in practice. If one has an ansatz for the dependence 
of the conformal moments on the complex conformal spin, one can employ for GDAs the partial 
wave expansion with integer conformal spin. 

In conclusion, we have introduced a representation that makes it easier to include GPDs and 
GDAs in phenomenological studies and offers new theoretical possibilities for the investigation of 
perturbative and non-perturbative aspects. 
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A Integrals 

In this appendix we collect several integrals, which appear in the evaluation of moments or the 
convolution of GPDs with the hard scattering amplitude. The partial waves Pj(x, rj), given in Eqs. 
(JHHJ), (J53|) . and (J55j), might be expressed by associated Legendre functions of the first and second 
kind, i.e., by a/(1 — x 2 )P J : I f 1 1 (x) and — x 2 )Q~ H 3 L 1 (x). The conformal moments Cj(x,rj), see Eq. 
(|HEJ), can be represented within the same basis, however, divided by the weight (1 — x 2 ) [75] . The 
integrals, presented in the following, can then be read off from diverse integral tables. Moreover, 
we give the relation between conformal and ordinary Mellin moments. 

The conformal moments of the partial waves for complex valued conformal spin read in the 
outer region 




(159) 
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where the normalization factor is given by 

Note that this integral (|159jl only converges for dtej > dtek. In the central region we have the 
following integral 

f^M- f'!ii + ')*W)M, (lei) 



'/ 



(A; + 1)0 + 2) 7T / j -k 

On the r.h.s. two terms appear in the brackets. The sum of the former one and the integral (|159|) 
is proportional to the identity S(j — k), however, the latter one gives an addendum. Fortunately, 
for integer value k = m = 0, 1, 2, • • • it will not contribute. Remarkably, if the support of Ck{x, 77) 
is extended to x < —rj by means of 



ri\2k+2 s [ n (k7t)T(l + k)r(3 + k) _ fe _ 3 /(fc + 3)/2,(fc + 4)/2 



T(3/2 + fc)r(5/2 + fc) z X V 5/2 + fc x 2 

and for pj(x,rj) by analytic continuation into the region x < 0, we have the following integral 

Jj*«S.*,nM*,n) = 1 1 gg£ 2) ^ffi? ■ »°*>«^- < 163 > 

Up to the sign it is equal to the sin(A:7r) proportional contribution on the r.h.s. in Eq. (jl61|) . Hence, 
the sum of integrals (|159|) . (|161j) and (|163j) might be understand as a limit 5ie k — > 9fte j that yields 
the identity, formally written as 

dxck(x,r))pj(x,r)) = — 2isin(irj)5(k — j) . (164) 

3 

We add that for non-negative integer values of the conformal spin the integral (|161|) establishes 
the orthogonality relation (|4*0|) for Gegenbauer polynomials. 

Let us also give here the Mellin moments of the partial waves pj (x, rj) for integer value n = 
0, 1, 2, • • -. For the integration in the outer region we have 

r, n , , 77"-Jsin(7T7) fl/2 + j/2,l+j/2,j/2-n/2 \ 

/ dxx Pj (x,ri) = — 3 F 2 ' ' 1 , $Kej > Ken, (165) 

J v (n-])ix V 5/2 + J, 1 +j/2-n/2 J 

while from the central region we find, up to the overall sign, the same expression 

, (j-n)ir 2 V 5/2+j,l+j/2-n/2 J 

Neglecting the sin(7rj) term, these expressions will contain single poles at j = {0,2, • ■ ■ ,n} and 
j = {1, 3, ■ ■ ■ , 71} for even and odd n, respectively. The values for the two lowest moments are: 

( , fA j 2 3 r(5/2 + j) sin(Trj) 
dxp j (x,r))= - . =t: , > , t , (167) 



V J jT(l/2)r(4 + i) 71 
= fgy ,. ,,™^L . (168 



rjj (j-l)(4 + j)r(l/2)r(3 + j) tt 
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Analogous as discussed above, the sum of both integrals (jl65j) and (j!66|) should be understood 
as a limit that results in a linear combination of "5-functions" , which are concentrated in j = 
{0, 1, • • • , n}. From the Mellin-Barnes representation of GPDs we explicitly find then the usual 
Mellin moments, expressed in terms of conformal ones: 



9 » m *_, l + (-l)"- t U!r(5/2 + z) 

dxX n q(x, 71, A 2 ) = V (i) , , V ; ; { , ■ ; — 771^7], A 2 ) . (169) 

.! yv ' ; ~^ 2y ^ 2z!r(l + n/2-z/2)r(5/2 + z/ 2 + ^/2) 
The inverse relation can be brought in the form 

/ a 2 \ fv\ 2i (-i)*n!r(3/2 + 7i-i) r 1 , „_ 2i , A2 , 

w n (?7, A 2 ) = g (-) l !( J_ 2 ^ (3/2 + n ; / 2 ^,,,A 2 ). (170) 

In the convolution of the hard-scattering amplitude with generalized parton distributions the 
following two integrals appear: 

00 dx . ^ [2\ 1+j r(5/2 + j) 

{ ^'«=y r(3/2)r(3+i) - (m) 



«-i-ie" v Vfy r(3/2)r(3+j) 

We note that the reduction to non-negative integer values of the conformal spin is simply done 
by the replacement j — > n. In this case only the central region contributes in these integrals, 
especially, the imaginary part will drop out 



B Gluonic sector 

Here we present the Mellin-Barnes integral for gluon GPDs, defined in Eq. (J2J). Let us first note 
that the index of Gegenbauer polynomials, appearing in the definition of conformal moments, is 
determined by group theory. To be more general, we consider the light-ray operator 

0(k 1 ,k 2 ) = (t){K 2 n)(j)(Kin) (173) 

that contains two quantum fields, which live on the light-cone n 2 = 0. For gluons the field is build 
by the field strength tensor. We assume that these fields have definite spin projection s on the 
light- cone, i.e., 

n^ v n v (j){Kn) = s <p(Kn) . (174) 

Here E^,, is the usual generator of Lorentz transformation, acting on a field 4>(x) at x — 0. 
Moreover, the canonical dimension of the field <f>(x) is denoted as I. The conformal spin of the 
field is defined as 

j = l(£ + s ). (175) 
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It characterizes the behavior of the field under collinear conformal transformation, which can be 
viewed as the projective transformations on a line: 

/ aK + b 

k — ► k = , ad — be = 1 , 

ck + a 

</>(kti) -> 0'M = (ck + d)- 2j <f> n ) . (176) 

V CK ~\- (X J 

In four dimensional space the quantum numbers are 

3 1 

£ = - , s = ±- for quark fields (177) 
£ = 2 , s = 0, ±1 for gluon field strength tensor (178) 

The conformal operators are obtained by the group theoretical decomposition of the light-ray 
operator ()173j) into irreducible representations. If both quantum fields have the same conformal 
spin, then the (local) conformal operators which have definite conformal spin are characterized by 
Gegenbauer polynomials with index v = 2j — 1/2: 

O nl = i l (d Kl + dj CI ( 9ki ~ M 0( Kl , k 2 ) , n<l, (179) 

where the conformal spin of these operators is 2j + n. For each given conformal spin, there appears 
an infinite tower of operators which are labelled by the quantum number I, related to their spin 
or if one likes to their canonical dimension. For leading twist operators the spin projection s must 
be maximal so that the twist of the fields t = £ — s is minimal. Consequently, to leading twist 
the quark fields have conformal spin j — 1 and gluon ones j = 3/2. The index of Gegenbauer 
polynomials is in the former and latter case v = 3/2 and v = 5/2. We remark that the gauge link 
factor connecting the fields along the light-cone does not change the construction via Eq. (|179|) . 
For gluonic GPDs (J2J) we define the conformal moments for n — 1, 2, ■ • ■ 

G„^^_n-iofA G„ /_\ _ r (5/2)r(n) 5/2 



c n (^v)=V n - LU c n ^-j with ^c n (x) = —^J^-^ C° n L\ (x) . (180) 

in such a way that in the forward limit the usual normalization of Mellin moments appear: 
lim^_>o c n{ x i v) = x n_1 • Moreover, as in the quark sector, they project on operators with confor- 
mal spin 2 + n. Compared to the Mellin moments of parton densities within our convention, one 
power in x seems to be missed, however, it is included in the definition of GPDs. For instance, in 
the vector case we have 

|m G /(i,^A 2 ) = ig(i) for x > , (181) 

where g(x) is the unpolarized gluon parton density. Analogous convention holds for the axial- 
vector case. The analytic continuation of the conformal spin j in the conformal moments is again 
done in terms of hypergeometric functions 
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The conformal partial waves are analogously constructed as described in Sect. 12.3.11 Namely, 
by including the weight (1 — x 2 /rf) 2 with a suitable normalization and from the requirement that 
the partial waves are vanishing at x — —77 and are continuous at the cross-over point x = rj: 

%(x, rj) = e(rj - \x\) V -t °r 3 (^j + 9(x - 77)77"* % (^j (183) 



where 



r(l/2)r(j) ( } 3 



x 



Gm ^ _ ^(Trj) r ^/2,(i + l)/2 



;i84) 
;i85) 



or 

Moreover, Bose symmetry implies definite symmetry of the gluonic GPDs (J2J) under the transfor- 
mation x — > —x, i.e., in the (axial-)vector case they are always (anti-) symmetric. This property is 
simply restored by forming symmetric or antisymmetric partial waves, see Eq. (|75jl . Corresponding 
to our normalization (j!81j) . we write 

1 rc+ioo (—1) 

G q(x, 77, A 2 ) = - / dj [ G p (x, 77) ± %(-x, t,)] mj (ri, A 2 ) , (186) 

2 J c _ ioo sm(TTj) 

where only the analytic continuation of conformal moments for even (odd) integer values of n is 
needed for (anti-)symmetric gluon GPDs. 

Let us comment on the properties of the gluonic conformal partial waves (|183|) . They are 
related to the quark ones by a derivation with respect to x: 

Pj(x,v) = ~-^Pi^ x ^) ■ ( 187 ) 

As a simple consequence, the first derivative of the gluonic partial wave is smooth at the cross-over 
point, while the second one has a jump as it is the case in the quark sector. Also at the point 
x = —rj the gluonic partial waves vanish as (x + rj) 2 rather than (x + 77) as for the quark ones. 
Consequently, the convolution with the following hard-scattering amplitudes, which appear in the 
electroproduction of transversely polarized vector mesons (70], exist 



(t + x)' \(J r(3/2)r(3 + J )' 1 ' 

iX %M = e-^(l) 1+ 'j^,: j) ,- (189) 



(£-x-iey- w r(3/2)r(3 +J -; 

It has been argued already in Ref. [7H] that this integrals without any further regularization exist, 
which is shown here from a more general point of view. We add that in the convolution of the 
gluonic GPDs with the hard-scattering amplitude at leading twist-two the integrals appear 

dx (i 9 o) 



* Z + x "» \a r(3/2)r(4 + J 



^- x - ie ^ Vey r(3/2)r(4 + J ) 



^.(x,0 = e-^- 1 M- — - ^— (191) 
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Finally, the convolution with conformal kernels is analogous done as outlined in Sect. |3] and 
yields in the Mellin-Barnes integral representation of GPDs to a multiplication of the conformal 
moments with its eigenvalues. Here we only mention that in correspondence with our normalization 
these eigenvalues follows from the Mellin transform of the forward limit, i.e., from the Mellin 
moments of the usual DGLAP kernels. 



C Mellin-Barnes representation of conformal kernels 

In the following we derive the Mellin-Barnes integral for a generic kernel (|57jl . which is conformal 
covariant and possess only non- negative eigenvalues k n . This is analogously done as for GPDs in 
Sect. 12 .HI Let us first remind on the support extension of the spectral representation (|87jl. The 
procedure for the kernel K(x, y) is well-known and arise from the representation 

K(x,y) = / dw + dw- 5(x — yw + — wJ) k(w + ,wJ) . (192) 

Jo J -l+w + 

Here k(w+, wJ) is the analog of the DD in the case of GPD, see Eqs. and (|21jl. Such a kernel 
appears in the convolution with a light-ray operator ()173)) and it is in the mathematical sense a 
distribution. The support in the whole (x, y)-plane can be read off from Eq. (|23|) and is written 
here as 

K(x, y) = 9(y - x)6{x + 1) [k(x, y) - 6{x - l)k(-x, -y)\ + ( X ~* ~ % \ , (193) 



y^-y 



where the distribution k(x, y) has the integral representation 



/ l+y 

k(x,y) = / dw + k(w + ,x — yw + ) . (194) 
Jo 

Analogous as in the case of a GPD, see discussion in Sec. 12.11 this integral is only uniquely 
defined in the central region, i.e., |x| < 1, while in the outer region 1 < |x| only the difference 
k(x, y) — k(—x, —y) enters. Here the ambiguities in the support extension of k(w+, x — yw + ) drop 
out. The support extension from the region \x\, \y\ < 1 to the whole support is unique |2], see 
analogous discussion as in the last paragraph of Sect. 12.21 

Now we are in the position to derive the Mellin-Barnes integral representation for the kernel 
K(x,y). To do so, we represent the series ()87|) in the region \y\ < 1 as integral in the complex 
plane that includes the positive real axis 

I r(°°) i 

K (%, y) = 7T 9 d J ■ r -\ PA x ) k 3 c Ay) ■ ( 195 ) 

2i / (0) sm(7rj) 

Here the functions Cj(x) = Cj(x, 1) and Pj(x) = Pj(x, 1) are defined in Eqs. (jSHj) and pHjl. respec- 
tively, and kj is the analytic continuation of the eigenvalues k n . They coincide with the Mellin 
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moments of the corresponding DGLAP kernel and might possesses a logarithmical growing for 
j — > oo. The integrand has simple poles at j — n — 0, 1, 2, • ■ • and so the residue theorem leads 
to the series (|87jl . Next, we deform the integration contour in such a way that it includes the 
imaginary axis and is closed by an arc with infinite radius, see Fig. |3] (a). It remains to show that 
this latter contribution vanishes. The behavior of the integrand for large j with |arg(j)| < ir/2 
can be estimated from the behavior of hypergeometric functions and is in our case given by 



-y g±j'{arccosh(— x— ie)±arccosh(y— ie)} 

■ f -s Pj(x + ie)cj{y ~ ie) ry— r . (196) 

sm(7rjj sm(7rj) 

Here the analytic continuation in the complex plane by the ie prescription determines the branch 
of the arccosh function, not uniquely defined for real valued x, y < 1. It is done in such a way 
that the estimate (J196)) is applicable for — 1 < y and x < 1 and |arg(j)| < tt/2. Obviously, on the 
infinite arc, both the denominator and numerator in Eq. (j!96|) will exponentially grow and the 
integrand will vanish as long as the condition |arccosh(— x — ie) ± arccosh(y — ie) J < n is satisfied. 
This is the case for — 1 < x < y < 1. Thus, for this region we arrive at the Mellin-Barnes 
representation for the Kernel K(x,y), i.e., for the function 

2 rc+ioo y 

2 Jc~ioo sm(7rj) 

Form the spectral representation ()192j) it follows that the convolution of the kernel with a holo- 
morphic test function t(x) yield a holomorphic function, depending on y. Thus, we might now 
employ analytic continuation to extend the representation (j!97j) into the region 1 < y. This will 
not alter the convergency properties of the integral, since arg (arccosh (y)) — for 1 < y. We 
remind that Pj{x) coincides for x < 1 with the holomorphic function Vj(x), defined in Eq. (JS2J). 
Vj(x) has a branch cut, starting at x — 1, along the positive real axis and we might define its 
value on the this cut by (Vj(x + ie) + Vj(x — ie))/2. Hence, within this procedure the function 
k(x,y) is uniquely continued into the region 1 < y for all values —1 < x, cf. with the prescription 

(USD- 

The Mellin-Barnes integral for the kernel K(x, y) follows now in an unique way from the result 
(j!97j) and its support, cf. Eq. (j!9Hj) . For 1 < x we can write K(x, y) as the difference 

K(x,y) = k(x,y) — k(—x, —y) for 1 < x < y . (198) 

We represent k(—x, —y) via the Mellin-Barnes integral (J197J) . where Pj{—x) is obtained by analytic 
continuation while Cj(—y) has now for 1 < y a branch cut and a single pole at y = 1 on the real 
axis. We have now to employ the principal value prescription, as above for Vj(x), and so its value 
on the cut is 

2 \ c i(-y + ie ) + c A-y - ie )\ = c °^M c j(y) + ^(^^djiy) , (199) 

where the function dj(y) is defined in Eq. 1)67)1 . Using the asymptotic behavior of hypergeometric 
functions for large j, it can be shown that in the region 1 < x < y the contribution Vj(—x)kjdj(y), 



59 



appearing in the integral (jl97j) . exponentially vanishes for j — > oo with |arg(j)| < tt/2. Hence, we 
can close the integration path so that the positive real axis is encircled. Since it is a holomorphic 
function, the corresponding integral vanishes. So we can drop this contribution and find that for 
1 < x < y the extension of the integral kernel, 



k- k- 
:Pj(x)cj{y) ~ 



sin(7rj) sin(7rj) 
of the Mellin-Barnes integral for K(x, y) 



1 1 

-Vj(x + ie) + -Vj(x - ie) - cos(7r jjVji-x) 



Cj (y), (200) 



1 rc+ico ^ 

2 Jc-ioc sm(Trj) 

The formula ()200|) defines the continuation of Pj(x) for 1 < x, which coincides with Eqs. fl53|) - (jo3|) 
for r\ — 1. The missing part of the kernel, i.e., for y < x, follows by the symmetry transformation 
x, y — > —x, —y from Eq. (J201j) . Hence, the integral kernel can be written as 



K{x,y) = ~ / dj^^r[9(y-x)p j (x)c j (y) + 0(x-y)p j (-x)c j (-y)}. (202) 

We excluded in our analyze here the line x = y. As long jkj is vanishing for j ; — ► oo, we can 
use analytic continuation to approach it. In all other cases a more advanced analyze is required. 
However, this is in fact not necessary here, since it is obvious that a constant or logarithmic 
behavior of kj for j ' ^ oo is associated with 5-functions and +-prescriptions for singularities at 
x = y. It is easy to check by forming the lowest moment that the Mellin-Barnes integral is correct 
in that case, too. 
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